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A GENERALIZED VORTEX LATTICE METHOD FOR SUBSONIC 
AND SUPERSONIC FLOW APPLICATIONS 

Luis R. Miranda, Robert D. Elliott, and William M. Baker 
Lockheed-California Company 


SUMMARY 


A vortex lattice method applicable to both subsonic and supersonic flow 
is described. It is shown that if the discrete vortex lattice is considered 
as an approximation to the surface— distributed vorticity , then the concept of 
the generalized principal part of an integral yields a residual term to the 
vorticity-induced velocity field. The proper incorporation of this term to 
the velocity field generated by the discrete vortex lines renders the present 
vortex lattice method valid for supersonic flow. Special technicjues for 
simulating nonzero thickness lifting surfaces and fusiform bodies with vortex 
lattice elements are included. Thickness effects of wing-like components are 
simulated by a double (biplanar) vortex lattice layer, and fusiform bodies 
are represented by a vortex grid arranged on a series of concentrical cylin- 
drical surfaces. The analysis of sideslip effects by the subject method is 
described. Numerical considerations peculiar to the application of these 
techniques are also discussed. A summary comparison of the results obtained 
by the method of this report with other theoretical and experimental results 
is presented. This method has been implemented in a digital computer code 
identified as VORLAX. A users manual for the VORLAX program is contained in 
Appendix A. A complete Fortran compilation and executed case are contained in 
Appendix B. Appendices C and D describe input conversion programs useful for 
transforming input between the VORLAX and NASA Wave Drag programs. 


INTRODUCTION 


The several versions or variations of the vortex lattice method that are 
presently available have proven to be very practical and versatile theoretical 
tools for the aerodynamic analysis and design of planar and nonplanar config- 
urations. The success of the method is due in great part to the relative 
simplicity of the numerical techniques involved, and to the high accuracy , 
within the limitations of the basic theory, of the results obtained. But 
most of the work on vortex lattice methods appears to have concentrated on 
subsonic flow application. The applicability of the basic techniques of 
vortex lattice theory to supersonic flow has been largely ignored. The method 
presented herein allows the direct extension of vortex lattice techniques to 
supersonic Mach numbers. The equations allowing this extended application are 
derived in the next section starting from the first order vector equations 
governing inviscid compressible flow. They are then applied to the particular 
case of a skewed-horseshoe vortex with special attention given to the super- 
sonic horseshoe. 


In the following theoreticsl discussion, the bssic srguinents involved in 
the simulation of thickness and volume effects by vortex lattice elements are 
presented. This particular modeling of the above effects represents an alter- 
native, with somewhat reduced computational requirements, to the method of 
quadrilateral vortex rings (refs. 1 and 2). The simulation of thickness and 
volime effects makes possible the computation of the surface pressure distri- 
bution on wing-body configurations. The fact that this can be done without 
having to resort to additional types of singularities, such as sources, 
results in a simpler digital computer code. ’ 
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THEORETICAL DISCUSSION 
The Basic Equations 


Weird has 
of an inviscid 
equations : 


shown, (ref. 3), that the small-perturbation, linearized flow 
compressible fluid is governed by the three first order vector 


V X V = U), 


V. w = Q, 


w = V (1) 


on the assumption that the vorticity oTand the source intensity_Q are toown 
functions of the point whose position vector is R. The vector v is 
perturbation velocity with orthogonal Cartesian components u, y, 

\ is a constant symmetrical tensor that for orthogonal Cartesian coordinates 
with the x-axis aligned with the freestream direction has the form 




2 

1 - M 0 

00 

0 1 

0 0 


0 

0 

1 


( 2 ) 


where is the freestream ^fech number. _If 0 = 
has the components w = 0^ui+vj+wk. This 
by Robinson (ref. 4 ), who called it the reduced 
denotes the total velocity vector, i.e. , u = (u» 
the fluid density, then it can be shown that for 


1-MoP, then the vector w 
vector was first introduced 
current velocity If u 
+ u) T + V J + w k, and p 
irrotational and horaentropic 


flow 


p u = Poo w + higher order terms 


( 3 ) 


where 
e.g. , uoo 
directly 


subscript «> indicates the value of the quantity at upstre^ infin y, 
= Uoo *i • Therefore, to a linear approximation, the vector w is 
related to the perturbation mass flux as follows: 


w = ( p u - 


u, 


o) /Pc 


( 4 ) 


The second equation of (l), i.e., the continuity condition, shows that for 
source-free flows (Q = O), w is a conserved quantity. 


Ward has integrated the three first order vector equations directly 
without having to resort to an auxiliary potential funct^n. He obtained 
two different solutions for v (r), depending on whether 0 is positive (su - 
sonic flow), or negative (supersonic flow). These two solutions can be com- 
bined formally into a single expression if the following convention is used. 


2 

K = 2 for 0 >0 

K = 1 for 0^ < 0 1/2 

Rg = Real part of |(x-xi)^ + 0^ [(y-y^) + J/ 
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Ilp = Finite part of integral as defined by Hadamard (refs. 5 
and 6 ) . 

The resulting solution for ttie per^rbation velocity v at the point 
whose position vector is i + y^ j + k, is given by 


’ (’* 1 ) ' - A 


2ttK 


~l~ n . V 


(R) ds 

% 


-j- I n X V ( R) [ 


R - R, 


dS 


^ -y- Q(R) V ^ dV -j~ -y- K i(R) dV ( 5 ) 

0 J R ^ 

V V 


This formula determines the value of v within the region V bounded by 
the surface S. The vector n is the unit outward (from the region V) normal 
to the surface S. Fiurthermore , it is understood that for supersonic flow 
only those parts of V ajid S lying within the domain of dependence (Mach 
forecone} of the point R^ are to be included in the integration. 

For source-free (Q = 0), irrotational C5=0) flow, equation ( 5 ) reduces 


-j- dS t^4-|Uy -1(8)1 ^dS( 6 ) 

s s ■'s 

This is a relation between v inside S and the values of n.w and n < 
on S, but- these two quantities cannot be specified independently on S. 

To determine the source-free, irrotational flow about an arbitrary body 
B by means of equation ( 6 ), assume that the surface S coincides with the 
wetted surface of the bo<^, with any trailing wake that it may have, and with 
a sphere of infinite radius enclosing the body and the whole flow field about 
It, namely, S = Sg -h S^, 

This surface S divides the space into two regions, Ve external to the 
body, and Vi internal to it. Applying equation (6) to both Ve and Vi, since 
the integrals over Sa, converge to zero, the following expression is obtained: 



vCRj^) = 


R— R 

-f N . (R)Vr ^ ^ jN X 6 v(R)t X ^ as (7) 


.here N = hi= -he 1= «>e uhlt normal to the body, or whe a_s the =at^. 
positive from the valu;s of the quantities on 

^hl ^oSe'spoiilng face f “rLrcfriatnStonTf s:r?:"1:Lity 

r;iT:Sfe\re%:rd1frf::e1nLri ilvea the contribution of a vorti- 
city distribution of surface density N x A V. 

If the boundary condition of zero mass flux through the surface 
Is applied to both external and internal flows 

N . p Ue = N . + Pco'^e^ = ° 

( S) 

N . p Ui = N . ( P„ + Poo ■”i^ ^ ° 

then the condition N . A ^ = 0 exists over Sg + S^, and the flow field is 


uniquely determined by 


Y (R) X -TT^ 




Where y (R) = ^ x a v is the surface vorticity density. 


Extension to Supersonic Flow 

In order to extend the application of thejortex lattice meth^ to 

process, the result being 

(ao) 


5 


where F= lim y • 6 


Y -*• oo 
6—0 

8 ia a dlMnaion norml to y, and di la the dlatance element along v. In the 
Classical vortex lattice method, applicable only to subsonic flo^^ the vorti- 
city distribution over the body and the wake, i.e., over the surface St. + Sr 

fields^*^^ ^ suitable arrangement of vortex filaments whose velocity 
fields are everywhere determined by equation (lO). This procedure is nJ 

supersonic flow. For this latter case, it is necessary 
donf in the foll^i^? approximation to it. This is 

o« surface Sg + Sw, which defines the body and its wake, is considered 

« being compoaed of a large number of diacrete fit area element r ^ef 

denaity y can be aasumed approximately eonatant 
then eolation (9) ean be approximated by the followlng^uatta: 


■ - h 


N 

Z 

J=1 


/ 


R-R, 


Xj ^ 


dS 


R' 


(11) 


e 


where N is the total number of discrete area elements T . When the point 
whose position vector is Ri is not part of Tj, the integral over this d^t, 
Crete area can be approximated by the mean value theorem as^followL- 


f 

Xr 


R-R, 




R- 


dS = Y, 


6j X 




/ 


R-R, 


R' 


di 


( 12 ) 


P 


T« ^ ’■J parallel to the average direction of y in T,, 5., 

Cj, and di is the arc length element along c/. Sis 
the velocity field induced by a discrete vorticity patch ^Tj can 
be approximated for points outside of Tj by some mean discretryortex'^line 
whose strength per unit length is yj But if the poinrif is 

of the^SuX integral in equation (ll) has an inherit sSgularity 

of the Cauchy type due to the fact that R = Rq at some point within r ll 

close to irblJt Lcat consider a point 

mite pert concept being eufflVlInt to LdSa^l^^S^.Sit^ 
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R-R, 


Y X ^ 
R^ 


dS = 


P 


lim 

e-.-o 


-L + lim -4- 

J e-»o J 


T- e 

R-R, 


= lim l( g) + y f> X 

G-»o 


JT «-«i 

J R^ 

n " 


di 


(13) 


The Oast integral in equation (13) repreaents “""J* 

vortex line contribution whose evaluation presents no difficulty. 
to determine the integration denoted by 

the coordinate system is centered at the point Ri, and tha ^ ^ - 

Ltermined by the discrete area t. Then, if Y denotes the modulus of Y, 


/ N _T y sin A - X cosA /p 


dx dy 


( 14 ) 


where A is the angle between the y-axis and the direction of the 

T!I t and b 2 = -s2 > 0 (supersonic flow). The components of the vector 

?orec<me from the point (o, o. t) with the T-plane, conaequently, if 
the integration with respect to x is performed first, 


l( e)=Y cosAj 
X 


f!f 


-B 


sf/ + ^ 

ty - X 


(y^ t 7W^ 


\ 

dx j dy 


ty -(1 +B)g 


where t = tan A , and Xi, \z are the wl^ues of y correspond! 
section of the line x=ty -(1 +B)g with the hyperbola 

2 

y » 


0 = G^(l-^2B)-2(l+B)Gty -(b‘ 

integration yields 


" t^) 


londi ng to 

= -B\/y2 + 


X = -BVjr + e 

then the finite part of the x- 


(15) 


to the inter- 
2. Let 


l( g) = Y cosA 


i 


ty (ty-(l+B)G) _ _1_ 

(y2 + 4T 


B 


! 


dy 


T 


V cos A 
■ b2 






2.2 

y + e 


_d 2 _ 

VX 


(16) 


Since e is a very small quantity, the variation of y in the interval ( Xj, X2) 
is going to be equally small, and, therefore, the quantity within brackets in 
the last integrand of eq^uation ( I6) can be replaced by a mean value and taken 
outside of the integral sign. The same is not true of the term since 

it will vary from 00 for y = Xq, go through finite values in the integration 
interval, and then again increase to cd for y = X 2 * With this in mind, and if 
y denotes a mean value of y, l( e) can be written as 


Ke) = 

B 


B^e-(l+B) cty - (B^-t^) y 


•>.2 2 

y + e 


/ 


i- 

v/^ 


(17) 


But \i, A 2 are the roots of ty-e = -BVy 
the polyncMlal denoted by Thus 


r~2 ^ 

V y + € , i . e . , they are the roots of 


x/c^(l+ 2 B)- 2 (l+B)ety - (B^-t^) y^ (l 8 ) 

Introducing this expression for into (I7), and taking the limit 
€-^o, the following value for l( r^) is obtained: 


^ 

I( o) = lim I( g) = - ^ — sj B^-t^ -j- ■ ■ ■ 

^ J x/Ui-yKy-X^) 


(19) 


The integral appearing in equation (19) can be easily evaluated by com- 
plex variable methods; its value is found to be 


f 


Xg 


dy 


\/( \j-y)(y-\2) 


= TT 


( 20 ) 
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The contribution 
•within T , induced by 
therefore given by 


of the inherent singularity to the velocity field, 
the vorticity patch T , and dented herein by t**, is 


« 


V 


lim I ( E ) 
2tt E o 


y cos A 
2 



( 21 ) 


This contribution is 
physical meaning when 
of the Mach lines. It is 
(11) and (12), that malies 


perpendicular to the plane of T, and it has only 
> t2 i.e., when the vortex lines are swept in front 
expreLion (2l), taken in conjunction with equations 
the vortex lattice method applicable to supersonic flow. 


The Skewed-Horseshoe Vortex 


Velocity fields due to complex vortex curves can ^^.^^neratedby the linear 

rsirvers“.‘a:; ‘ ° 


u 


+ 


C!o ^2 

2itK 



dy 

I (x-x^)2 (y-y^)^ + ) 


3/2 


V = 


C!o p2 

2vK 


w = 


+ 


r 

2vK 




(x-x ) dy - (y-y^) ^ 

I {K-xf » (y-y„)^ + 


( 22 ) 
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where P represents the circulation per 
length, and the integrations are to be 
satisfies the conditions 


unit length of discrete vortex line 
carried out along that part of C which 


(x - X ) +p {(y - y )^ + z ^) > 0 


and 


X - X < 0 
o 


if 


> 1 . 


For the transverse leg of the horseshoe vortex, the coordinate x appearing 

indlrat^r"^ expressed as a function of y, i.e., x = ty, and the 

indicated integrations carried out between the limits y = -s and y = +s 

figure 2* By defining the following auxiliary variables * 


x_ = X + ts 
1 o 


7l = yo + s 


= X - ty 
o o 


X - X - ts 

2 o 


= 7o " ^ 


pie resulting formulas giving the velocity components induced by the skewed 
transverse rectilinear vortex filament can be written as follows: 


u = + 


Fz . 

2ttK ’ ^2 

x + 


(t^+|3^) 


tXi + ^ 


tx2 + ^ y. 


_Vxi^ + 


V = - 


2ttK 




* 0 


* P Yr 


yj\ ~P 


w = - 


2vK 


+ (t^+|3^) z ^ 
^ O 


^^1 ^1 


tx„ + /J y. 




ft.l/Lw ®^ressions, the coordinates (x„, yo, Zo) °f the receiving 

ItlLf tl with respect to the midpoint of the rectilinear vortex 

coordinate system coinciding with the plane de- 
fined by the x— axis and the vortex itself. 
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( 23 ) 


i 



The case of a rectilinear vortex segment parallel to the x-axis (t = « ) as 
of special importance since the trailing legs of a horseshoe vortex are gen- 
erally asstamed to be parallel to the x-axis. Since dy = 0, equations (22) 
■become 


u = 0 


V = - 


2ttK 


t 


dx 


|(x-x^)^ + + z/)} 


3/2 


(2U) 




dx 


2ttK 




If the vortex segment extends from x - x^ to x - Xj., the above integra- 
tion yields 


u = 0 


p2 

I O 

2ttK 


X - X. 

o 1 










(25) 


W = + 


2ttK 


X - X, 
O . 1 


X - X- 
o f 


^(Xq-x.)^ + Uy^-yf * z/) ^(x^-x^)^ + ((y^-y 


) . z/) 


y 




For a conventional horseshoe vortex whose trailing legs stretch to down- 
stream infinity, equations (25) would give the contribution of the port leg with 
the following substitutions 


x. = « if M^<1 

1 


X. 

1 


Xo - ((y^+s)^ + z/) 


i/ >1 


X 


f 


-ts 


y = -S 
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. Likewise, the contribution from the starboard trailing leg can be computed 
by introducing the following values into equations (25) 


= +ts 

= “ if < 1 

Xf = x^ - if M„>l 

y = +s 


Combining these results with equations (23), the formulas defining the 
now field induced by a discrete vortex consisting of a skewed segment and two 
railing legs parallel to the x-axis (the skewed-horseshoe vortex) are obtained. 
Keeping in mind Hadamard's finite part concept, and after introducine the fol- 
lowing notation 

F = '^’‘2 * 

(26) 

°1- 'r , * C = 1; M. >1: C = 0) 

V«i *(! (yi 


I 2 2 2 T" ^ ^ = LJ >1: C = 0 ) 

V-2 


then, the horseshoe vortex induced field formulas can be expressed as follows: 
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“ ^o' 


= + 




(Fi - F^) 


'' ‘‘o’ ^o’ 


= + 


(F^ - F^) t 


2ttK o j x*^ + (t^ + /3^) 7.^ y^ + 


2 2 2 2 


(27 


^2 ■" "o 


^ ^o’ 


X (F^ - F^) 


G, - 


y/*z/ ^ 


2 o 


The finite part concept determines the value of the constant C appearing in 
the definition of and G^- 

A notable simplification of equations (27) occurs for supersonic flow 
when the receiving point (x^, y^, z ) is in the plane of the horseshoe, namely, 
Zq = 0- First, the values of the axialwash and sidewash, u and v, vanish 
identically; secondly, the upwash expression becomes 



Equation (28) is applicable to both subsonic and supersonic flow in its present 
format. But for the supersonic flow case, the fact that the constant C of the 
G functions becomes null due to the finite part concept allows further simpli- 
fication of the upwash equation. Introducing the corresponding values of the 
F and G functions, the expanded version of equation (28) is 


w (x , y , 0) 
o o 


2ttK 





L 

tXi + 

tx^ + P^Y2 

■■ 

/ 2 ^ ,2 2 

i 2 „2 2 

+ P 

\ 


*i/yi 


Xj./y2 


y/ ^Xj^ * y 


2 2 
2 


(29 
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Since x = 


= X. 


- ty^ = x-L - tyj - X2 - ty2, the rearrangement of equation (28) 


/ 


in factors of |^(xj_ - ty^) v'xi 
finally reduces the upvash formula to 


'Z d ~ ^ \ y 

^2 yi2j-l and |^(x2 - ty2) ^ 


w (x^, y^ 0) = 


2itK 




+ P y 


2 2 
2 


(30) 


When the field point (x^, y^, 0) is within the distributed vorticity patch 
whichis approximated by the discrete horseshoe vortex, e.g. , the control point 
associated with the horseshoe, the upwash given by equation (30) has to be 
complemented by the distribution due to the generalized principal part of the 
upwash integral, as givenby equation (21). If 6x is the distance, measured 
in the X direction, occupied by the distributed vorticity Y which has been 
liunped into the discrete transverse vortex leg of circulation R the relation- 
ship between the Y of equation (21) and the Rof equations (27) and (30) is 


r = Y cos A 6x 


(31) 


Modeling of Lifting Surfaces with Thickness 


The method of quadrilateral vortex rings placed on the actual body sur- 
face (ref. 1) provides a way of computing the surface pressure distribution 
of arbitral^ bodies using discrete vortex lines only. Numerical difficulties 
may occur when the above method is applied to the analysis of airfoils with 
sharp traili^ edps due to the close proximity of two vortex surfaces of 

parallel direction. An alternative approach, requiring somewhat less 
counter storap and easier to handle numerically, consists in using a double 

horseshoe vortices to model a lifting surface 
with thictaess, as shown schematically in figure 3. This constitutes an 

location of the singularities, similar in nature to 
classical lifting surface theory approximation of a cambered sheet. 

swept horseshoe vortices, and their boundary condition control 
points, corre^onding to a given surface, upper or lower, are located in a 
same P^ne. The upper and lower surface lattice planes are separated by a 
gap which represents the chordwise average of the airfoil thickness distri- 
bution. The results are not too sensitive to the magnitude of this gap; any 
TO ue between one half to the full maximum chordwise thickness of the airfoil 
has been found to be adequate, the preferred value being two thirds of the 
mximum thickness. Furthermore, the gap can vary in the direction normal to 
the x-axis to anow for spanwise thickness taper. On the other hand, the 
chordwise distribution, or spacing, of the transverse elements of the horse- 
shoe vortices have a significant influence on the accuracy of the computed 


li* 



surface pressure distribution. For greater accuracy, for a given chordwise 
number of horseshoe vortices, the transverse legs have to be longitudinally 
spaced according to the cosine distribution law 



( 32 ) 


where xj - Xq represents the distance from the leading edge to the midpoint 
of the swept leg of the Jth horseshoe vortex, c is the length of the local 
chord running through the midpoints of a given chordwise strip, and N is the 
number of horseshoe vortices per strip. The chordwise control point location 
corresponding to this distribution of vortex elements is given by 


The control points are located along the centerline, or midpoint line, of the 
chordwise strip (fig. H). Lan has shown (ref. 7) that the chordwise 'cosine' 
collocation of the lattice elements, defined by equations (32) and (33), 
greatly improve the accuracy of the con 5 >utation of the effects due to lift. 

His results are directly extendable to the computation of surface pressure 
distributions of wings with thickness by the biplanar lattice scheme pre- 
sented herein. 

The small perturbation boundary condition 

V . n' = -u . n (3^) 

00 

is applied at the control points. In equation (3^)» n = il + mj + nit, and 
n^ = mj + lik, where i, m, and n are the direction cosines of the normal to 
the actual airfoil surface. Equation (3^) implies that llu| « | mv + nw | • 

The vise of the small perturbation boundary condition is consistent with the 
present biplanar approach to the simulation of thick wings. 


1 - cos 


H) 


(33) 


Modeling of Fusiform Bodies 


The modeling of fusiform bodies with horseshoe vortices requires a 
special concentrical vortex lattice if the simulation of the volvime displace- 
ment effects, and the computation of the surface pressure distribution, are 
to be carried out. To define this lattice, it is necessary to consider first 
an auxiliary body, identical in cross-sectional shape and longitudinal area 
distribution to the actual body, with a straight baricentric line, i.e., 
without camber. The cross-sectional shape of this auxiliary body is then 
approximated by a polygon whose sides determine the transverse legs of the 
horseshoe vortices. The vertices of the polygon and the axis of the auxiliary 
body (which by definition is rectilinear (zero camber) and internal to all 
possible cross sections of the body) define a set of radial planes in which 
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the bound trailing legs of the horseshoe vortices lie parallel to the axis 
(fig. 5). As the body cross section changes shape along its length, the 
corresponding polygon is allowed to change accordingly, but with the constraint 
that the polygonal vertices must always lie in the same set of radial planes. 
The axial spacing of the cross-sectional planes that determine the transverse 
vortex elements, or polygonal rings, follows the cosine law of equation ('^ 2 ). 
The boundary condition control points are located on the auxiliary body sur- 
face, and in the bisector radial planes, with their longitudinal spacing given 
by equation (33). ^ & o 

The boundaiy condition to be satisfied at these control points is the 
zero mass flux equation. 


w.n=-u.n (35) 

OC 

where all the components of the scalar product w . n = u + mv + nw 
are to be retained. Thus, equation (.35) is a higher order condition than 
equation (3^0. The use of this higher order boundary condition, within the 
framework of a linearized theory, is not mathematically consistent. There- 
fore, it can only be justified by its results rather than by a strict mathe- 
matical derivation. In the present treatment of fusiform bodies, it has been 
found that the use of higher order , or exact boundary conditions is a re- 
quisite for the accurate determination of the surface pressure distribution. 


The fact that the vector w, instead of v, appears in the left hand member 
of equation (35) requires some elaboration. First, it should be pointed out 
that for small perturbations w , n^ v . h'. Furthermore, for incompressible 
flow (B = 1), the vector w is identical to the perturbation velocity v. Con- 
sequently, the boundary condition equation (3^^) is consistent with the con- 
tinuity equation, V. w = 0, to a first order for compressible flow, and to any 
onder for incompressible flow. But when a higher order boundary con- 
dition is applied in compressible flow to a linearized solution, it should be 
remembered that this solution satisfies the conservation of "w, not of v, i.e., 
V. w = 0. Thus, the higher order boundary condition should involve the 
reduced current velocity, or perturbation mass flux, vector w, as in equation 
( 35 ), rather than the perturbation velocity vector 3F. 

The body camber, which was eliminated in the definition of the auxiliary 
body, is taken into account in the computation of the direction cosines i,m, 
and n, which are implicit in equation ( 35 ). Therefore, the effect of camber 
IS represented in the boundary condition but ignored in the spatial placement 
of the horseshoe elements. This scheme will give a fair approximation to 

cambered fusiform bodies provided that the amount of body camber is not too 
large* 
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Computation of Sideslip Effects 


The aerodynamics of an isolated wing in sideslip can be analyzed by two 
different approaches depending on the coordinate system chosen. In one approach, 
the coordinate system consists of wind axes, the longitudinal axis being aligned 
with the free-stream velocity vector, figure 6. This formulation of the problem 
is known as the skewed-wing approach, and a first order solution obtained within 
this framework will give the dominant effects of sideslip, even for the case of 
zero dihedral. The other approach, also shown schematically in figure 6 and 
known as the skewed free-stream approach, is based on a body-axis formulation 
of the problem and the corresponding first order solution, though it may be 
adequate for large dihedral, will fail to produce the significant effects of 
sideslip for low or zero dihedral. To compute the sideslip effects correctly 
within the framework of a skewed free-stream formulation, it is necessary to 
solve partial differential equations containing second order terms of the per- 
turbation velocities. This implies a much more involved computational procedure 
than that required for the solution of the first order perturbation equations (l) 
On the other hand, the application of the skewed-wing approach to anything more 
complex than an isolated wing in sideslip, such as might be the case with a 
configuration with wing, fuselage, and nacelles, becomes geometrically very 
complicated. 

The approach adopted herein is a combination of the two approaches men- 
tioned above, formulated with the objective of obtaining reasonably accurate 
sideslip effects using only a first order perturbation solution but without all 
of the geometrical complications inherent in the skewed-wing approach. Basically 
it is assumed that the vortex lattice representing the configuration and its 
vortex wake consists of both bound and free elements or legs; the vortex fila- 
ments that model rigid surfaces are considered bound, and those that constitute 
the wake are the so-called free elements, figure 7* The bound portion of the 
lattice, containing both transverse and trailing, or chordwise, segments, is 
invariant in a body axis system, the chordwise legs being parallel to the x-axis. 
The free legs of the lattice are not actually force free, rather they are as- 
sumed to extend to downstream infinity parallel to a predetermined direction 
Ts proportional to the angles of attack and sideslip. If the propor- 
tionality factors are unity, then the free portion of the lattice would be 
invariant in wind axes . 

After the circulation strengths of the above lattice geometry are solved 
for under the appropriate boundary conditions, the pressure coefficient dis- 
tribution is computed in accordance with the higher order expression 

c = - (0^ u + V„ v) (36) 

q^c 

where Uoo and Voo are the components along the x and y body axes of the free-^ 
stream velocity vector of modulus ; the corresponding perturbation velocity 
components are denoted by u and v, as usual. The use of equation (36) instead 
of the linear approximation 
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Uoo U 


(37) 


is required for the correct computation, within the present framework, of the 
rolling moment due to sideslip of a planar wing. This is due to the fact that 
the hound trailing legs, being defined in body axes, are not lined up with the 
free-stream flow and therefore, according to the theorem of Kutta, they con- 
tribute to the normal force. This contribution is represented by the second 
order term in equation (36), namely, v. Even though this contribution is 
of second order, it must be included in the computation of the differential 
rolling moment due to sideslip, since this quantity itself is of the second 
order for a planar, or nearly planar, wing. 


THE GENERALIZED VORTEX LATTICE METHOD 
Description of Computational Method 

The four items discussed in the preceding section, i.e., the inclusion of 
the vorti city-induced residual term w* for supersonic flow, the biplanar scheme 
for representing thickness effects, the use of a vortex grid of concentrical 
polygonal cylinders for the simulation of fusiform bodies, and the special 
lattice geometry consisting of both bound and free elements for the analysis 
of sideslip effects, have been implemented in a computational procedure herein 
known as the generalized vortex lattice (GVL) method. This method, outlined 
in what follows, has been codified in a Fortran IV computer program (VORLAX). 

The basic element of the method is the swept horseshoe vortex whose trail- 
ing legs has both bound and free segments. The latter segments may trail to 
downstream infinity in any arbitrary, but predetermined, direction whereas the 
bound trailing legs are laid out on the proper cylindrical control surfaces in 
a direction which is parallel to the x body axis. Figure 8 illustrates sche- 
matically the representation of a wing-body configuration within the context 
of the present method. In this illustration, the streamwise arrangement of 
the lattice follows the cosine distribution law, equation (32), but both chord- 
wise and spanwise distributions of vortex lines can be independently specified 
to be either of the cosine or of the equal spacing. To each horseshoe vortex 
there corresponds a control point which is placed midway between the bound 
trailing legs of the horseshoe; the longitudinal location of the control point 
is determined by equation (33) if the cosine chordwise distribution has been 
chosen, otherwise it is located halfway between the transverse legs, as 
required by quarter-chord/three-quarter-chord rule. 

The direction of floatation of the wake vortex filaments is defined by 
the two angles and shown in figure 8, the former being proportional to 
the angle of attack, and the latter being proportional to the sideslip angle. 
The proportionality constants are part of the program input, the recommended 
values being 0 for the sideslip constant, and 0 or 0.5 for the angle of attack 
constant. 
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The velocity field induced by the elementary horseshoe vortex is given by 
(27) wLn above constants are both sero, and by somevhat more compll- 
cS e^reslions which take into account the hinks in the ‘-“I-* 
either one or both of the wake floatation parameters are nonzero. Though not 
preLnted here, these expressions can be easily derived through the application 

of equations (23)* 

When the influence induced by a horseshoe vortex upon its own control 
point is being evaluated, the contribution from the generalized ’ 

at; eiven bv equation (2l), is included in the normalwash if the free-stre^ is 
sUrsonic! . for the supersonic case, the simplified downwash 

fomula, equation (30), is used whenever the receiving point is in the same 
plane of the inducing horseshoe. 

The horseshoe vortex velocity field is used to generate the coefficients 
of a system of linear equations relating the unknown vortex strengths to the 
:LropriIt: bLdary condition at the control points Ihls ™ 

sSved by either a Gauss-Seidel iterative procedure, known as controlle 
successive over-relaxation (ref. 8), or by a vector orthogonalization tech- 
nique! ^.e.. Purcell's vector method (ref. 9). If the inverse process is 
desired, i.e., synthesis or design instead of analysis, the linear sys em 
eruaiions is ^sed to compute the slope distribution 

to achieve a specified surface loading; this involves a straightforward matrix 
m^tlpiCation process. Mired cases, i.e., design and analysis, are easily 
handled by proper grouping of the boundary condition equations. 

The pressure coefficients are computed in terms of the perturbation 
velocity components, the computation being carried out according o 
one of three possible ways, as follows: 

1 If the surface under consideration is assumed wetted by the flow 
on both sides (zero thickness panel) and the configuration si e- 
slip angle is zero, then a net loading coefficient is computed 
based on the local value of the spanwise vorticity, namely, 

ACp = 2 y cos A; 

2. When the configuration sideslip anjle is not 

coefficient of a zero- thickness surface is calculated through the 
use of the higher order expression (36); and 

3 When surface pressure coefficients are computed, i.e., thepanel 
under consideration is assumed wetted by the flow on one side 
only, the isentropic flow relationship giving the pressure 
coefficient in terms of free-stream Mach number and local-to free- 
stream velocity ratio is resorted to. 

The force and moment coefficients are calculated by nmerical Integra^ 

of the pressure coefficient distribution with <iue VORLAX 

edge forces. If cosine chordwise lattice spacing is specified by the VORU 
program user, the computation of the leading edge suction of zero thickness 
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panels is carried out according to Lan's procedure (ref. 7 ), whose application 
to supersonic flow is made possible by the generalized vortex-induced velocity 
field formulas presented in this report. If equal chordwise lattice spacing 
is specified, the contribution of the leading edge suction singularity to the 
forces and moments is calculated by the technique indicated by Hancock in 
reference 10; this approach is not nearly as accurate as Lan's, the magnitude 
of the leading edge suction being significantly underestimated. 

The VORLAX computer program has the capability of analyzing symmetrical 
and asymmetrical cases as well as configurations in steady state angular rota- 
tion about any or all of three axes, parallel to the coordinate axes, going 
through the input moment reference center. Steady state angular rotation cases 
are treated by the subterfuge of assuming a nonuniform onset flow, this onset 
flow being defined by the values of the angular rates and distance of the field 
point to the rotation center. 

^Ground proximity effects are analyzed by the method of im.ages , i.e., the 
configuration is mirrored about the ground plane; the flow around the airplane 
and its image then contains a stream surface which coincides with the ground 
plane due to the symmetry of the arrangement. In this modeling of a configu- 
ration in ground proximity, it is assumed that the trailing vorticity wake 
floats to downstream infinity parallel to the plane of the ground. 


Numerical Considerations 

At supersonic Mach numbers, the velocity induced by a discrete horseshoe 
vortex becomes very large in the very close proximity of the envelope of Mach 
cones generated by the transverse leg of the horseshoe. At the characteristic 
envelope surface itself, the induced velocity correctly vanishes, due to the 
finite part concept. This singular behavior of the velocity field occurs only 
for field points off the plane of the horseshoe. For the planar case, the 
velocity field is well behaved in the vicinity of the characteristic surface. • 
A simple procedure to treat this numerical singularity consists of defining 
the characteristic surfaces by the equation 


(x-x^)^ = C |(y-y3^)^ + (z-Zj^)^j- (38) 

where C is a numerical constant whose value is greater than, but close to, 1. 

It has been found that this procedure yields satisfactory results, and that 
these results are quite insensitive to reasonable variations of the paran- 
eter C. 

Another numerical, prohlem, peculiar to the supersonic horseshoe vortex, 
exists in the planar case (field point in the plane of the horseshoe) when the 
field point is close to a transverse vortex leg swept exactly parallel to the 
Mach lines (sonic vortex), while the vortex lines immediately in front of and 
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behind this sonic vortex are subsonic and supersonic, respectively. Ih 
problem can be handled by replacing the boundan' condition ecuati.'m for 
sonic vortex with the averaging equation 


uch 


-1 +1 


where Yt* is the circulation strength of the critical horseshoe vortex, and 
and Yi*+i respective circulation values for the fore-and-aft 

adjacent subsonic and supersonic vortices. 


Tlie axialwash induced velocity component (u) is needed for the computation 
of the surface pressure distribution, and for the formulation ol the boundary 
condition for fusiform bodies. When the field point is not too close to the 
generating vorticity element, the axialwash is adequately' described by the ^ 
conventional discrete horseshoe vortex representation. But if this point is 
in the close vicinity of the generating element, as may occur in the biplanar 
and in the concentrical cylindrical lattices of the present method, tiie error 
in the computation of the axialwash due to the discretization of the vorticity 
becomes unacceptable. This problem is solved by resorting to a vortex-splitting 
technique, similar to the one presented in reference il. Briefly, this tech- 
nique consists of computing the axialwash induced by the transverse leg of a 
horseshoe as the summation of several transverse legs longitudinally redis- 
tributed, according to an interdigitation scheme, over the region that contains 
the vorticity represented by the single discrete vortejc. This is done only if 
the point at which the axialwash value is required lies within a given near 
field region surrounding the original discrete vortex. 


COMPARISON WITH OTHER THEORIES AND EXPERIMENTAL RESULTS 


Conical flow theory provides a body of exact results, within the context 
of linearized supersonic flow, for some simple three-dimensional configurations. 
These exact results can be used as bench mark cases to evaluate the accuracy of 
numerical techniques. This has been done rather extensively for the GVL method, 
and very good agreement between it and conical flow theory has been obsei’ved in 
the computed aerodynamic load distribution and all force and moment coefficients 
Only some typical comparisons are presented in this report, figures 9 through 12 

Finally, the capability of computing surface pressure distributions by the 
1110 'thod of this paper is illustrated in iigures 13 and lU. 
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CONCLUDING REMARKS 


The present vortex lattice method, in the form of a computer program, 
has the capability to calculate the aerodynamic load distribution at subsonic 
and supersonic Mach numbers for arbitrary nonplanar configurations. It has 
been found to be a very useful preliminary design tool, particularly when 
aircraft configurations whose mission requirements involve both subsonic and 
supersonic flight are considered. It is also capable of the inverse process, 
nmely, the computation of the surface warp required to achieve a given load* 
distribution. Correlation with experimental data and with results from other 
theories shows a good agreement not only in the overall force and moment 
coefficients due to lift, but also in the distribution of the load coefficients. 

The schemes shown for the simulation of thickness and volume effects, 
which allow the computation of surface pressure distribution by using only 
vortex lattice singularities, appear adequate for most practical purposes, 
though experience in this respect is somewhat limited. 

The treatment of sideslip cases by the present method does not require 
higher order solutions, as is necessary for the skewed free-stream approach, 
and it is not as geometrically complicated as the skewed-wing formulation. 

Yet the analysis of complex configurations in sideslip still requires care 
and caution due to the numerical anomalies that may result from the inter- 
action among aircraft components, such as a horizontal' tail or a body, and 
the "free" trailing legs of the horseshoe vortices . 

Additional capabilities that can be added to the present computer code, 
and that would enhance the value of the method as a preliminary design tool 
include the following: 

• Incorporation of an optimization algorithm based either on Lagrange 
multipliers or on a gradient method, to design the surface warp 
for minimum drag under specified constraints. 

• Application of the technique of reference 11, or of some other 
adequate technique, for the simulation of jet exhaust effects, 
with particular attention to its extension to supersonic flow. 

• Introduction of a design procedure for the calculation of the 
geometry required to achieve a given surface pressure distribution, 
i.e., synthesis of both camber and thickness. The biplanar vortex* 
lattice simulation of a thick lifting surface is well suited for 
the development of such a design procedure when combined with an 
iterative scheme. 
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Figure 1.- Definition of integration regions for the 
computation of principal part* 
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Figure 2.- Sweptback 
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Figure 3.- Modeling of thick wing with horseshoe vortices. 



Figure U.- Vortex lattice collocation 






BOUNDARY CONDITION: 
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Figiire 5*- Modeling of fusiform body with horseshoe vortices. 
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HORSESHOE FREE LEGS 
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Figure 8.- Generalized vortex lattice model of wing-body configuration- 






Figure 9.- Theoretical comparison of arrow wing lift slope and aerodynamic center location 



Figure 10.- Theoretical comparison of arrow wing drag-due-to-lift factor 




Figure 11.- Theoretical comparison of chordwise loading for delta wing. 


SWEEP=59“i ASPECT RATIO = 1.92 MACH NUMBER=v/2 
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Figure 12. Theoretical comparison of chordwise loading for sveptback rectangular ving 
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Figure 13.- Comparison with experimental pressure distribution 
on wing-body model at Mach = 0.5. 













Figure Ih.- Theoretical comparison of pressure distribution on 

ellipsoids at zero angle of attack in incompressible 
flow 


APPENDIX A 


USER'S MANUAL 
FOR 

A GENERALIZED VORTEX LATTICE METHOD 
FOR SUBSONIC AND SUPERSONIC FLOW APPLICATIONS 



THE VORLAX COMPUTER PROGRAM 


A computer program has been developed for the aerodynamic analysis ^d 
design of arbitrary aircraft configurations in subsonic and supersonic flow. 
This computer program, herein identified as VORLAX, has been codifie in 
FORTRAN IV for use on the CDC 6600, the IBM 360, and the jy 

computer systems. A complete compilation and executed case in CDC F^RT:^ I 
is contained in Appendix B. Two auxiliary interface programs, 'treated in 
Appendices C and d! provide for input data transformation from NASA Wave Drag 
format (Reference 12) to VORLAX format (Appendix C) and for transformation 
from VORLAX to NASA Wave Drag format (Appendix D) . 

The VORLAX program is based on a generalized vortex lattice (GVL) method 
Which extends the applicability of vortex lattice .techniques to 
range of problems than has heretofore been considered. In this program, the 
configuration is represented by a three-dimensional, generally nonplanar, 
vortex lattice; the basic element of the lattice is the skewed horseshoe 

vortex whose induced velocity formulas have ^dLSe^Mplanar 

and supersonic flow. Thickness effects can be simulated by a ^ 

or sandwich) lattice arrangement. Fusiform bodies can be modelle y 
centric cylindrical lattice of polygonal cross sections . The computational 
capabilities of the program include the following: 


• Surface pressure or net load coefficient distribution. 

• Aerodynamic force and moment coefficients. 

• Surface warp (camber and twist) design in order to achieve 
input pressure distribution, 

• Longitudinal/lateral stability derivatives • 

• Ground and wall (wind tunnel interference) effects. 

• Flow field survey. 

• Symmetrical/asymmetrical configurations and/or flight conditions 

The limitations of the VORLAX program are characteristic of methods 
based on inviscid linearized potential flow theory, as follows; 


• Attached flow. 

• Small perturbation flow, 

• Flow entirely subsonic or supersonic (no mixed transonic flow) 

• Straight Mach lines . 

• Rigid vortex wake. 
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PRACTICAL INPUT INSTRUCTIONS 


In defining a configuration for the program input, a master frame of 
reference X- Y- Z- is assumed. The X-Z plane is the centerline plane, the 
Z axis directed upward, and the X-axis pointing in the downstream direction; 
the Y axis points to starboard. The origin of the system can be any con- 
venient point in the X-Z plane. In general, the configuration can be made 
up of symmetrical and asymmetrical components, and in defining the sym- 
metrical components only the starboard elements need be specified. 

The configuration to be input is divided into a set of major panels; 
up to 20 of these panels can be input, symmetrical components being counted 
only once. For instance, a wing with straight leading and trailing edges, 
and with linear lofting between the root and tip, constitutes a major panel. 
Complex planforms , and nonlinear changes in twist and airfoil sections are 
described by defining more than one panel for a given wing. The computer 
program will then subdivide each major panel into a number of smaller 
elem.entary panels spaced chordwise and spanwise, i.e., a finer mesh lattice 
is generated internally. The chordwise and spanwise spacing is specified 
by the user, two options being available: (l) the semicircle or cosine 

distribution so well known in airfoil and wing theory, and (2) the equally 
spaced distribution. 

Up to 2000 elementary panels can be used in the definition of a given 
configuration. Any consistent system of length and area dimensions can be 
used in the specification of the configuration length, but it is recommended 
that the system of units used be such that the largest length dimension does 
not require more than three digits to the left of the decimal point. Other- 
wise, significant digits may be lost in the output printout. 

Wing thickness effects can be taken into account within the context of 
control surface theory . By control surface theory one means that the exact 
linearized theory is used to evaluate induced velocities along a given mean 
surface, known as the control surface, and these values enter into the com- 
putation of the boundaiy conditions which are satisfied at this control surfaci 
rather than at the actual boundary surface. The control surface equivalent 
of a typical two-dimensional airfoil is illustrated in figure A-1. The 
assumption inherent in control surface theory is that the induced velocities 
vary very little in the vicinity of the surface. Experience has shown that 
for the majority of practical configurations the loss in accuracy is neg- 
ligible, and is more than compensated for by the increase in computational 
efficiency. Any wing-like component with thickness is then represented by 
a double set of panels, one for the upper surface and one for the lower 
surface as sho-.”- schematically in figure A-2. 

All the swept horseshoe vortices, and their boundary condition control 
points, corresponding to a given surface, upper or lower, are located in a 
same plane. The upper and lower surface lattice planes are separated by 
a gap which represents the chordwise average of the airfoil thickness 




hi 


Figure A-1. - Control surface nomenclature 
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Fi^Jiure A-2. - Modeling of thick ving vith horseshoe vortices 


distribution. The results are not too sensitive to the magnitude of "this 
gap- any value between one half to the full maximum chordwise thickness of 
the airfoil has been found to be adequate, the preferred value being two 
thirds of the maximum thickness. Furthermore, the gap can vary in the 
direction normal to the x-axis to allow for spanwise thickness taper. On 
the other hand, the chordwise distribution, or spacing, of the transverse 
elements of the horseshoe vortices have a significant influence on the ac 
curacy of the computed surface pressure distribution. For greater accuracy, 
for a given chordwise number of horseshoe vortices, the transverse legs 
have to be longitudinally spaced according to the cosine distribution law 



where xY - represents the distance from the leading edge to the midpoint 
of the swept leg of the Jth horseshoe vortex, c is the length of the local 
chord running through the midpoints of a given chordwise strip, and N is the 
number of horseshoe vortices per strip. The chordwise control point location 
corresponding to this distribution of vortex elements is given by 


- *0 - I [i - i )] 


The control points are located along the centerline, or midpoint line, of 
the chordwise strip (figure A-3). 

The modeling of fusiform bodies with horseshoe vortices requires a 
special concentric vortex lattice if the simulation of the volume displace 
ment effects, and the computation of the surface pressure distribution, are 
to be carried out. To define this lattice, it is necessary to consider 
first an auxiliaiy body, identical in cross-sectional shape and longitudinal 
area distribution to the actual body, with a straight bancentric line, i.e., 
without camber. The cross-sectional shape of this auxiliary bodj is t 
approximated by a polygon whose sides determine the transverse ° 

horseshoe vortices. The vertices of the polygon and the axis ol tne aux 
iliary body, which by definition is rectilinear (zero camber) and internal 
to all possible cross sections of the body, define a set of radial Pl^nes 
in which the bound trailing legs of the horseshoe vortices lie parallel to 
the axis (figure A-k) . As the body cross section changes shape along 
length, the corresponding Wygon is allowed to change accordingly, but with 
the constraint that the polygonal vertices must always lie in the same set 
of radial planes. As in the case of the biplanar representation of thinness 
effects, cosine ax'ial spacing should be used for the analysis of ;^sifora 
bodies. The effect of body camber is taken into account by independently 
specifying the camber of the baricentric axis of the body. 




Figure A-3. - Vortex lattice collocation 
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AXIAL DISTRIBUTION OF B.C. CONTROL POINTS: 


INPUT CARD IMAGE DESCRIPTION 


The cards that constitute the input data deck are described in the 
following paragraphs. All format specifications are given in FORTRAN IV. 

While some 19 cards are described, this does not mean that input for a 
case consists of I9 cards. Rather, these should be thought of as card types. 
Furthermore, not all card types will be included in a given case, those to 
be included or deleted being a function of some of the input values as shown 


CARD 

(1): 

TITLE 

CARD 

(2) : 

IS0LV 


LAX 


LAY 


In colTimns 1 through 80 write any alpha- 
numeric identification heading. 

Method to be used for solving the system 
of linear equations relating the boundary 
conditions to the vorticity strength, in 
integer format, in coliamn 2. IS0LV = 0: 
Gauss-Seidel relaxation with accelerated 
convergence (under- or over-relaxation); 

IS0LV = 1: Purcell's vector orthogonali- 

zation method. 

Chordwise, or streamwise, spacing of vor- 
tices, integer quantity punched in column 12, 
LAX == 0: vortices are collocated at the 

percent chord (X/c) values determined by 
the cosine law X/C - 0.5 ( 1-cos (( 2K-l) 
i"/2N)), where K varies between 1 and N, N 
being the number of chordwise vortices; 

~ 1" vortices are collocated according 
to the equally- spaced quarter-chord law 
x/c = (4K-3)/(Un). The cosine law is 
recommended for greater accuracy for a given 
number of vortices. 

Spanwise, or lateral, spacing of vortices, 
integer quantity punched in column 22. 

LAY = 0; vortices are spaced at intervals 
(elementary vortex span) given by the cosine 
distribution law Ab = 

(jTf/M))/2, where Ab is 

span, bp is the panel span, and J varies 
between 1 and M, M being the spanwise number 
of vortices in a given panel; LAY = 1: 
vortices are equally spaced along the panel 
span, i.e., Ab = bpy^pj. The cosine spacing 
is recommended for enhanced accuracy, but 


bp( cos( ( J-l) tt/m) -cos 
the vortex element 


h6 








REXPAR 


HAG 


FL0ATX 


FL0ATY 


for most cases the difference in the results 
between cosine spanwise spacing and even 
spanwise spacing is negligible. 

Over-relaxation parameter, in FIO.O format 
starting in column 31. This parameter is 
intended to accelerate the Gauss-Seidel 
relaxation process, and/or make it convergen 
when it might otherwise diverge. Blank, 
or zero, input, implies that the program 
will compute internally the optimum over- 
relaxation value. If a positive quantity 
between 0.01 and 0.99 is input, this becomes 
the value of the over-relaxation parameter 
that the program will use, the optimum value 
being overridden. If IS0LV = 1, this 
parameter is not used, and therefore, not a 
required input quantity. 

Height above ground of the moment reference 
center, in FIO.O format starting in column 
Ul. If it is punched equal to zero, or 
left blank, the height above the ground is 
infinity, i.e., no ground effect. If a 
quantity different than zero is input, then 
the ground effect will be computed by the 
method of images, the height being given by 
the input value, in consistent units. 

Longitudinal vortex wake flotation factor, 
in FIO.O format, starting in column 51. If 
zero, or blank, then the trailing vortex 
legs being shed from the corresponding 
trailing edges, extend to infinity parallel 
to the X-Y plane. If a value different froir 
zero is input, then the trailing vortex legs 
shed from the trailing edges form an angle 
»^=FL0ATX- ALPHA with the X-Y plane, where 
ALPHA is the freestream angle of attack. 

(See figure A-6) 

Lateral vortex wake flotation factor, in 
FIO.O format, starting in column 6l. If 
zero or blank, then the trailing vortex 
legs being shed from the corresponding 
trailing edges extend to infinity parallel 
to the X-Z plane. If a value different 
from zero is input, then the vortex legs 
shed from the trailing edges form an angle 


48 



HORSESHOE FREE LEGS 



Figure A-6. Generalized vortex lattice inodel of ving body configuration 




= FL0ATY-BETA with the X-Z plan, where 
BETA is the freestream angle of slideslip 
(see figure A-6). 

ITRMAX Maximum number of iterations allowed for 

the Gauss-Seidel relaxation method, in 13 
format right-adjusted to column 80. If no 
value is input, the code will make ITRMAX 
= 99 by default. If IS0LV =1, i.e., the 
vector orthogonal! zat ion solution is re- 
sorted to, then ITRMAX is not a required 
input . 

CARD (3): NMACH Number of Mach niimbers to be analyzed, in 

12 format; i.e., integer value of NMACH in 
column 2. NMACH <7* 

MACH Mach numbers in FIO.O format starting in 

column 11. 

CARD (i+): NALPHA Number of angles of attack in 12 format; 

i.e., integer value of NALPHA in coluinn 2. 
NALPHA < 7. 

ALPHA Angles of attack in degrees in FIO.O format, 

starting in column 11. 

CARD(5): LATRAL Asymmetric flight or configuration flag. 

0 in column 2 = symmetric flight and sym- 
metric configuration about the X-Z plane. 

1 in column 2 = asymmetric flight and/or 
asyramet r i c con f i gur at i on , 

PSI Sideslip angle in degrees in FIO.O format 

starting in coliunn 11. + = wind coming from 
left side of nose. Input 0. or blank when 
LATRAL is 0. Used to obtain static stability 
derivatives such as Cyp, etc. 

PITCHQ Pitch rate in degrees /second in FIO.O format 

starting in column 21. + = nose up pitch. 

Used to obtain dynamic stability derivatives 
such as Cm^, Cl^^, etc. LATRAL may be 0 
when PITCHQ is nonzero. 

R0LLQ Roll rate in degrees /second in FIO.O format 

starting in column 31. + = left roll. 

Input 0. or blank when LATRAL is 0. Used 
to obtain dynamic stability derivatives 
such as C^ , ^n-py etc. 
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YAWQ 


VINF 


CARD ( 6 ) : NPAR 


SREF 


GEAR 


XBAR 


ZBAR 


Yaw rate in degrees/second in FIO.O format 
starting in column ^tl. + = left yaw or air- 
stream component washing from left to right 
across nose of airplane. Input 0. or blank 
when LATRAL is 0. Used to obtain dynamic 
stability derivatives such as * C > etc. 

r •'r 


Reference free stream velocity in FIO.O 
format starting in colvunn 51- If value 
is input, VINF is automatically set equal 
to 1.0 by the program. This parameter is 
only used when any of the angular rates is 
different from zero. It enters in the 
computation of the equivalent flow angle. 

For instance, if VINF = WSPAN/2 (wing semi- 
span) and ROLLQ = 5*73, then pb/2V = 0.1, 
and the rolling moment coefficient printed 
out by the program will be exactly one-tenth 
the value of the stability derivative . 
Likewise, if VINF = CBAR/2 (half the meah 
aerodynamic chord) aiid PITCHQ = 5.73> then 
the difference between the output pitching 
moment coefficient and the pitching moment 
coefficient for the case PITCHQ = 0. will 
be equal to one-tenth of the C^ derivative. 

q 

Number of major panels that will define the 
configuration, in 12 format; i.e., integer 
value of NPAN in columns 1-2 right adjusted 
to 2. NPAN <20. 


Reference area for force and moment coeffi- 
cients, in FIO.O format starting in coltimn 
11 . 


Pitching moment coefficient reference length, 
in FIO.O format starting in column 21. 

Usually mean aerodynamic chord length. 

Abscissa of moment reference point, in 
FIO.O format starting in column 31. 
X-coordinate in master frame of reference. 

Ordinate of moment reference point, in 
FIO.O format starting in column Ul. 
Z-coordinate in master frame of reference. 
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WSPAN 


CARD (7): XI 


Y1 


Z1 


C0RD1 


CARD ( 8 ) : X2 


Y2 


Z2 


C0RD2 


Total wing span in units consistent with 
SREF and CBAR, in FIO.O format starting in 
column 51. If left blank, a value of 2.0 
will be assum.ed by the program. 

X or longitudinal coordinate of the leading 
edge of one side of a major panel. Usually 
taken as the most inboard side in the case 
of wings. Input in FIO.O format starting 
in column 1. 

Y or lateral coordinate of leading edge of 
first side of a major panel. Input in 
FIO.O format starting in column 11. 

Z or vertical coordinate of leading edge of 
first side of a major panel. Input in FIO.O 
format starting in column 21. 

Chord length of first side of major panel 
measured from XI, Y1 , Z1 above in the 
positive direction of, and parallel to, the 
X axis. 

X or longitudinal coordinate of the leading 
edge of the second side of the major panel 
described on card (7). Usually taken as the 
most outboard side in the case of wings. In 
the case of a closed curved panel, e.g. , a 
cylindrical segment representative of a 
nacelle, X2 would be identical to XI. Input 
in FIO.O format starting in column 1. 

Y or lateral coordinate of leading edge of 
second side of the major panel described on 
Card ( 7 ). Input in FIO.O format starting 
i^ column 11. 

Z or vertical coordinate of leading edge of 
second side of the major panel described on 
Card ( 7 ). Input in FIO.O format starting 
in column 21. 

Chord length of second side of major panel 
measured from X2, Y2, Z2 above in the posi- 
tive direction of, and parallel to, the 
X axis. 


Note; Columns hl-80 of cards ( 7 ) and (8) are not read. Thus, any informatior 
useful for identification purposes may be written there. 
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Note that the side edges 1 and 2 of a major panel, which have just heen^ 
described by the input in cards (7) and (8), define the direction oi ^the posi- 
tive normal for that panel. This is determined by the feet-to-head direction 
for an observer standing on the panel looking upstream (down the negative 
X-axis) and with panel edge 1 to his left, and panel edge 2 to his right. In 
this analogy it is assumed that gravity is not a factor in that an observer 
could be standing on the bottom of a wing, for example, equally as easily os 
he might stand on top. In the case of curved panels, such as a nacelle com- 
ponent, the direction of the positive normal loses its meaning for the panel 
as a whole, but it is still applicable to each one of the longitudinal, or 
chordwise, strips that make up the major panel. 

CARD (9) NV0R Number of spanwise elements or vortices that 

will be used to represent the p^^nel , in FIO.O 
format starting in column 1. KV0R £100.0. 

pjjQY Number of chordwise vortices that will be 

used to represent the panel, in FIO.O format 
starting in column 11. RliCV £50. Tne 
program, using NV0R and RI.CV, will subdivide 
the panel under consideration into a grid of 
NV0R X RNCV swept horseshoe vortices collo- 
cated in accordance with the values of the 
LAX and LAY param.eters. Note that cliordwise 
and lateral distributions are independent, 
e.g., a cosine chordwise spacing (LAX^O) is 
compatible with equal spanwise distribution 
(LAY = l), and vice versa. The corresponding 
control points at which the boundary condi- 
*^Lons are satisfied are collocated according 
to the law (X/C)control = U-cos (Ktt/K)), 

K varying between 1 and N, where N is the 
number of chordwise vortices, if LAX = 0. If 
LAX =1, then the control points are placed 
at (X/C)control = ( ^K-1 ) / ( ^0 , namely, accor- 
ding to the equally spaced three-quarter 
chord distribution. The spanwise location 
of the control points is always at the center- 
line of the elementary swept horseshoe 
* vortices. 


To determine the surface slope at each ele 
ment control point, the program uses straight 
line element lofting between the two longi- 
tudinal, or butt line, edges of the major 
panel . 
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In addition to the above limitations to the 
values of NV0R and MCV, if more than one 
major panel is used in the description of 
the configuration, the following should be 
observed : 

NPAN 

If LATRAL=0, S NV0R x RNCV <2000 
1 

NPAN 

If LATRAL=1, S NV0R x RNCV x I QUANT <2000 
1 

IQUANT being defined further down. 

SPC Leading edge suction multiplier, in FIO.O 

format starting in col^unn 21, 0. = no 

suction. 1. = 100 percent leading edge 
suction. Nonzero values are recommended for 
all panels whose leading edges are vretted by 
the airstream. The program has the capability 
of computing the effects of free leading edge 
vorticity (leading edge vortex flows) by a 
localized application of the Polhamus analogy. 
This computation is triggered by inputting 
the SPC parameter as a negative quantity. 

When this is done, the sectional leading edge 
suction vector will be rotated normal to the 
camber surface at the leading edge, instead 
of the corresponding attached-flow tangential 
orientation, and the forces and moments will 
be computed using the rotated suction vector. 

PDL Planar/curved panel flag, in FIO.O format 

starting in column 31. 0. = planar major 

panel is to be described (including warped 
planar.) PDL = 999- (or >360.) = a curved 
major panel is to be described. 

CARD (lO): PHI Polar coordinate angle of radius vector when 

defining the subpaneling of a curved major 
panel. Omit this card when PDL - 0 . PHI is 
the angle measured from the horizontal in a 
plane parallel to the Y-Z plane. PHI = 0 
coincides with a line parallel to and in the 
positive direction of the Y-axis. Positive 
values of PHI are measured counterclockwise 
when viewed from the rear of the aircraft. 
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R0 


CARD (ll): AIKCl 


AINC2 


ITS 


PHI is input in FI 0.0 format starting in 
column 1. This and the subsequent input RO 
constitute a polar coordinate pair. The 
number of pairs to be input = NV0R + 1. Four 
pairs per card may be input in FIO.O format 
starting in column 1. As many cards as 
necessary are used. While the location of 
the origin from which the polar angle, PHI, 
is arbitrary the first PHI, R0 polar coordi- 
nate pair must coincide with the Yl, Z1 rec- 
tangular coordinates input on Card ( 7 ) • 
Likewise, the last input polar pair must 
coincide with the Y2, Z2 of Card (8). 

Radius vector from arbitrary origin when 
defining the subpaneling of a curved major 
panel. Input in FIO.O format. Each R0 is 
part of a PHI, R0 polar coordinate pair. 

Tangent of the angle subtended by major panel 
root chord] ine, or first edge (described in 
Card (7)), and the positive X-axis, in FIO.O 
format starting in column 1. Sign convention 
is determined by observing the edge 1 chord- 
line and the X-axis from edge 2. The edge 1 
chord is then rotated counterclockwise until 
it is parallel to the X-axis. If the angle 
rotated through is less than 90 degrees then 
the angle, and consequently its tangent, are 
considered positive. If it is greater than 
90 degrees, then AINCl is negative. 

Tangent of the angle subtended by major panel 
tip chordline, or second edge (described in 
Card (8)), in FIO.O format starting in column 
11. Sign convention is determined by 
observing edge 2 and the X-axis , looking in 
the direction from edge 2 toward edge 1. The 
edge 2 chord is then rotated counterclockwise 
until it is parallel to the X-axis. If the 
angle rotated through is less than 90 degrees, 
then the sign is positive; otherwise it is 
negative. 

Surface flag input as a two place integer in 
columns 21 and 22, right— ad justed to column 
22. ITS = 0 or blank indicates that the 
panel is considered as a lifting surface of 
zero thickness, i.e., both its upper and 
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NAP 


I QUANT 


ISYNT 


lower surface are wetted ly the external flow* 
ITS = 01 means that only the panel upper 
surface is wetted hy the real external flow. 
ITS = -1 indicates that only the panel lower 
surface is wetted by the real external flow. 

A double panel setup can then be used to 
represent a wing-like component with non-zero 
thickness, as previously illustrated in 
Figure 2. Notice that the XI, Yl, Zl, X2, 

Y2, and Z2 values to be input correspond to 
the control surface plane, and not to the 
actual chordal plane. The results are not 
critically sensitive to the separation between 
the upper surface panel (ITS = 01 ) and the 
lower surface panel (ITS = -1 ) , a separation 
of two thirds the thickness ratio of the air- 
foil being a good average value to use. 

Niomber of percent chords or stations along 
the chord (C0RD1 and C0RD2) at which the 
camber, or surface, ordinates are to be 
input. Input as a two-place integer in 
columns 31 and 32, right-adjusted to 
column 32. Maximum value of NAP is 50. A 
NAP = 0, 1, or 2 will be interpreted as a 
flat wing and no subsequent camber cards will 
be expected. If ISYNT, on this same card, is 
to be input as 1, i.e., a design case, then 
NAP should be 0 or blank. 

Symmetry flag with respect to X-Z plane input 
as an integer in column k2. IQUANT = 0 or 2 
indicates there is a mirror image of the 
panel on the opposite side of the X-Z plane. 
IQUANT = 1 indicates the panel is unique to 
the side for which it is being input. 

Design/anal^^sis flag input as an integer in 
coliomn 52. ISYNT = 0 or blank indicates that 
the panel has been defined geometrically and 
only analysis is to take place. ISYNT = 1 
indicates that the panel camber is to be 
designed by the program to support a speci- 
fied pressure distribution. If NAP on this 
same card was input >2, then ISYNT should 
be zero or blank. 
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NPP Nonplanar parameter, input as on integer in 

column 62. NPP = 0 indicates that all the 
vortex filaments representing a given surface 
lie in the cylindrical surface whose directrix 
is the leading edge of the panel, and whose 
generatrices are all parallel to the X-axis, 
NPP = 1 denotes that the transverse vortex, 
filaments are located on the actual "body 
surface, hut the hound trailing legs are par- 
allel to the x-axis. This parameter affects 
the definition of ZC^ and ZC2 on cards 16 
and 18. 

CARD (12): Cl Pressure coefficients along the first, or 

root, edge of major panel defined on Card (7). 
If ISYNT = 0, then this card is omitted . 

Units are dimensionless , Ap/q. Format is 
8F10.0 starting in column 1 with as many 
cards as necessary to input RNCV values of Cl. 

The desired values of the aerodynamic loading 
are defined at the chordwise location of the 
vortex lines. Thus, if LAX =0, the cor- 
responding X/C points follow the cosine 
distribution (l-cos ( (2K-1)it/ 2N) )/2; if 
LAX = 1, then the definition points are 
located hy the law (^K-3)/^N. In the above 
expressions K ranges between 1 and N, N 
being equal to RNCV, the chordwise number of 
vortices . 

CARD (13): C2 Pressure coefficients along the second, or 

tip edge of major panel defined on card (8). 

If ISYNT = 0, then this card is omitted . 

Format is 8 FIO.O starting in column 1 with 
as many cards as required to input RNCV 
values of Cl. Linear interpolation between 
corresponding values of Cl and C2 is used to 
obtain AC at intermediate spanwise values 
for the su?panels. 

CARD (li+): XAF Chord percent values at which camber, or sur- 

face ordinates will be supplied for the major 
panel, in 8F10.0 format starting in column 1 
using as many cards as necessary to define 
NAP values of XAF. If NAP is 0, 1, or 2, 
then a flat uncambered surface is implied and 
this card is omitted. These chord percents 
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need not be equally spaced, but the same set 
applies to both the root and tip chords, or 
edges of the panel. Second order Lagrange or 
sliding parabola interpolation is used to 
interpolate between input XAF points to 
obtain the surface slope value at the control 
point of each subpanel. 

For a panel representing the surface of a 
non-zero thickness airfoil or body, ITS 5^ 0, 
the Lagrange interpolation is modified to a 
fractional power (l/2) Lagrange method in the 
neighborhood of the leading edge. This allows 
the precise definition of the surface slopes 
at the control points for a blunt leading 
edge. 

CARD (15): RLEl Leading edge radius, in percent chord, of 

airfoil section at the first, or root, edge 
of panel, in FIO.O format between columns 1 
and 10* This card exists only if ITS ^ 0 , 

PPL <360»0 and NAP >2, i.e., an airfoil with 
non-zero thickness is being simulated* Other- 
wise, it must be omitted* 

CARD (l6): ZCl Camber ordinates or surface ordinates of root 

chord of the major panel described on Card (7) 
in units of percent chord, in 8F10.0 format, 
using as many cards as necessary to input NAP 
values, each corresponding to an XAF of the 
Card (lU) series. Omit this card if NAP is 
0> 1^ or 2 * If a lifting surface, such as a 
wing, is being simulated by a zero thickness 
panel then the ordinates of the wing camber 
line should be input here* If the wing thick- 
ness is being simulated by the sandwich or 
biplanar, method, i.e., a separate panel for 
upper and lower surfaces, then the surface 
ordinates for upper or lower should be input 
here* If a curved panel is being simulated 
(i.e., PDL >360 .) and NPP = 0 then ordinates 
of the panel streamwise edge should be input. 
These might represent the mean line of a 
shaped cowl for a flow through nacelle, for 
example. In this simulation all shed vor- 
tices will lie in the same cylindrical sur- 
face determined by the leading edge of the 
curved panel and the X direction. If NPP = 1 
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and PPL > 36O.Q, the ZCl array represents the 
camber of the body axis. 

CARD (17): RLE2 Leading edge radius, in percent chord, of 

airfoil section at the second, or tip, edge 
of panel, in FIO.O format between columns 1 
and 10. This card exists only if ITS 0, 

PDL < 360.0 and NAP > 2, otherwise it must be 
omitted. 

CARD (18): ZC2 Camber ordinates, or surface ordinates of the 

second, or tip chord of the major panel, or 
area ratios of the major panel in 8F10.0 
format, using as many cards as necessary to 
input NAP values, each corresponding to ar 
XAF of the Card (l4) series. Omit this card 
if NAP is 0, 1, or 2. Linear spanwise inter- 
polation is used to obtain intermediate values. 

If a curved panel is being simulated 
(PDL > 360. ) and NPP = 1 then area ratios in 
percent are expected here. These should 
represent the ratio of the cross sectional 
area of the closed polygonal surface being 
simulated at the XAF station under considera- 
tion divided by the area of the reference 
polygon times 100. Note that ZC2 is entered 
in percent where a value of 100 represents a 
section exactly the size of the reference 
polygon. Values greater or less than 100 
are permitted down to and including 0. The 
reference polygon is that input via the PHI-RO 
pairs on Card 10. In this simulation it is 
presumed that all stations along the panel 
have the same shape as the reference polygon, 
and the transverse vortices are located on 
the actual body surface of the curved panel. 

This concludes the input for the first major panel of the configuration. 

If there is more than one panel, then start over with Card (?) and work down 
to this point. Panels may be input in any sequence. 

After the last panel is described, continue with Card (19) . 

CARD (19): NXS Number of X— stations that will define the 

spatial flow field survey grid. NXS = 00 
means no survey desired. Maximum value of 
NXS is 20. Input as a two-digit integer in 
columns 1 and 2, right-adjusted to column 2. 
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NYS 


NZS 

CARD (20): XS 

CARD (21): YN0T 

DELTAY 

ZN0T 

DELTAZ 


Number of Y-stations that vill define the 
butt line values of the survey grid. NYS = 00 
for no survey. NYS and NZS (following) may 
be any positive integer subject to NXS x NYS x 
NZS <2000. 

Input as a two digit integer in columns 11 
and 12. 

Number of Z-stations that vill define the 
water line values of the survey grid. 

NZS = 00 for no survey. Input as a two digit 
integer in columns 21 and 22. 

X station values for the spatial flow field 
grid. Omit this card if NXS=0 . Input in 
8F10.0 format starting in column 1, using as 
many cards as necessary to define NXS values. 

Beginning of grid in the butt line direction. 
Input in FIO.O format starting in column 1. 
Omit this card if NXS = 0. 

Y-spacing of the grid. There will be NYS 
butt line planes equally spaced a distance 
DELTAY apart. Input in FIO.O format starting 
in column 11. 

Beginning of grid in the water line direction. 
Input in FIO.O format starting in column 21. 

Z-spacing of the grid. There vill be NZS 
water line planes equally spaced a distance 
DELTAZ apart. Input in FIO.O format starting 
in column 31. 


This ends the input description for a single case. 

Consecutive data sets or cases can be submitted at the same time. The 
program will always identify the presence of a new set by the corresponding 
title card (Card (l)). 
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PROGRAM OUTPUT 


The program output is processed hy a standard 132 characters-per-line 
printer. The output from each configuration is preceded by a printout of the 
input data cards. This printout is not an exact image of the input deck; 
rather, it is the version of the deck as the code sees it, namely, the default 
value of an input parameter is printed if there is a corresponding blank in^ 
the input card. Also, data within a format field are lined up for clarity in 
identification, even though in the input deck such data may be arbitrarily 
located within its field. The input deck data is followed by a list of the 
major geometric parameters used by the program and generated from the input 
data deck. Next, the component and total force and moment coefficients are 
printed out for a given flow condition (Mach number, angle of attack, angle 
of sideslip, and rotational velocities). These are followed by a tabulation 
of the location of all the vortex elements, the pressure coefficients, the 
circulation strengths, and other ancillary information. If a flow field 
survey about the configuration has been requested, then the flow parameters 
(velocity components, flow angles, Mach number, and pressure ratios/ at a 
series of field grid points will be listed. If other flow conditions have 
been analyzed, the same type of output will follow for each one of them, 
starting with the listing of the component and total force and moment coeffi- 
cients. If other configurations have been input, then the output will continue 
with the listing of the corresponding input data deck, and so on. Rather than 
describing the output format in detail, a glossary of the output terms, 
arranged in sequential order of appearance, and a sample computer output. 

Table A-1, are presented. 

I Numbering index for major panel identification. For cases where 

LATRAL = 1, the I-number preceded by a double asterisk in the 
PANEL GEOMETRY list denotes that the panel is the mirror image 
(about the X-Z axis) of the panel with the same I-number but 
without asterisks. 

XAPEX =XI (see input terminology). 

YAPEX =Y1 (see input terminology). 

ZAPEX =Z1 (see input terminology). 

PDC =PDC (see input terminology). 

LESWP Sweep of the panel leading edge, in degrees. Positive for 

sweepback, negative for sweep forward. 

CSTART =C0RD1 (see input terminology). 

TAPER Panel taper ratio, C0RD2/C0RD1. 
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PSPAN Panel span . 

NV0R =NV0R (see input terminology). 

RNCV =RNCV (see input terminology). 

SPC =SPC (see input terminology). 

SURF Panel surface area. 

CN Panel normal force coefficient, referenced to its own siirface 

area, 

CL Panel lift coefficient, wind axes, referenced to its own 

surface area. 

CY Panel lateral force coefficient, wind axes, referenced to its 

own surface area. 

CD Panel drag coefficient, wind axes, referenced to its own 

surface area. 

CT Panel leading edge thrust coefficient, referenced to its own 

surface area. 

CS Panel leading edge suction coefficient referenced to its own 

surface area. 

CM Panel pitching moment about moment reference center divided 

by (freestreajn dynamic pressure X SURF), wind axes. 

CRM Panel rolling moment about moment reference center divided by 

(freestream dynamic pressuire X SURF), wind axes. 

CYM Panel yawing moment about moment reference center divided by 

(freestream dynamic pressure X SURF), wind axes. 

SREF =SREF (see input terminology). 

WSPAN =WSPAN (see input terminology). 

CBAR =CBAR (see input terminology). 

CLT0T Total (summation over all panels) lift coefficient referenced 

to SREF. 

CDT0T Total pressure drag coefficient, referenced to SREF. 

CYT0T Total lateral force coefficient, wind axes, referenced to SREF. 
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CMT0T Total pitching moment coefficient about moment reference center, 

wind axes, referenced to SREF and GEAR. 

CRT0T Total rolling moment coefficient about moment reference center, 

wind axes, referenced to SREF and WSPAN. 

CNT0T Total yawing moment coefficient about moment reference center , 

wind axes, referenced to SREF and WSPAN. 

E Oswald's efficiency factor. 

S Perimetral, or spanwise, index of vortex element. 

C Chordwise, or streamwise, index of vortex element, 1 denotes 

leading edge element, and C value equal to RNCV corresponds to 
the last, or trailing edge element. 

X/C Percent chord location of boiuid. vortex line. 

X, Y, Z Coordinates of horseshoe vortex element centroid (center point 

of bound vortex line). 

CH0RD Local chord length. 

SL0PE Surface slope at boundary control point. 

ITS Flag which indicates type of panel surface (see input terminol- 

ogy). A zero value means that the panel is considered as a zero 
thickness lifting surface. A positive unit value (l) denotes 
that the panel is the upper surface of an airfoil-like element. 

A negative unit value (-1) corresponds to the lower surface. 

In the case of body-like components, 1 denotes the external, or 
wetted, surface. For a flow-through nacelle arrangement, 1 
stands for the external surface, and -1 for the internal surface. 

DCP Local loading coefficient (ACp = Cpg - if panel ITS = 0. 

If ITS ^ 0, then DCP is the local pressure coefficient (Cp). 

CNC Sectional normal force coefficient times local chord. 

CN Sectional normal force coefficient. 

DL Local dihedral, in degrees. 

CMT Sectional pitching moment coefficient about local quarter chord. 

GAMMA Vortex element circulation strength, divided by freestream velocity. 

ZC/C If the design option is being invoked, the resulting surface warp is 

printed out instead of GAMMA. This surface warp is expressed in 
fraction of the local chord, and it includes both camber and twist. 

CTC Sectional thrust coefficient times local chord. 
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CDC 

ITRMAX 

EPS 

ITER 

BIG 


Sectional pressure drag coefficient times local chord. 

Maximum allowable number of relaxation cycles if Gauss-Seidel 
solution is used (see input terminology). 

Tolerance, or minimum, iteration change for relaxation process. 
Actual number of relaxation cycles. 

Actual iteration change for relaxation process. Relaxation will 
stop when BIG < EPS, or when ITER = ITRMAX. 


if 8- flow field survey about the configuration has been requested 

(NXS > 0 in the input deck, see input description), then the following param- 
eters will also be output as part of this survey: 

X, Y, Z Coordinates of the survey grid points (not to be confused with 
the X, Y, Z coordinates of the vortex centroids previously 
described), referenced to the configuration master coordinate 
•frame. The X, Y, Z coordinates are determined by the XS, YNOT, 
DELTAY , ZNOT , DELTAZ input values (see input terminology). 

U Dimensionless velocity ( freestream velocity at infinity assumed 

to be unity) along the X-direction (body axes). 

V Dimensionless velocity along the Y-direction (body axes). 

W Dimensionless velocity along the Z-direction (body axes). 

EPSL0N Upwash angle in degrees, this angle is measured with respect 

to the X-axis in a plane parallel to the X-Z plane. 


SIGMA 


Sidewash angle in degrees. This angle is measured with respect 
to the X-axis in a plane parallel to the X-Y plane. 


CP 

MLOC 


Local field pressure coefficient (CP 
Local Mach number. 


^^static"^inf ^/'^inf ^ ' 


P/PT0T Local static-to-total pressure ratio for isentropic flow. 

P/PINF Local to-freestream static pressure ratio for isentropic flow. 
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TABLE A-1. SAMPLE OUTPUT 

IirUT DECK 


iftiditc Cl^UtSC VEHICLE CCLl607-3?/-?4 
0 I 0*10 


t 

0.500 

2.550 


1 

5.000 



• 

0.0 

0.0 

O 

. 

o 

11 

9676.80 

102.2280 

191.0400 



>444 


0.0 

0.0 

4.5762 

352.3999 

0*0 

0.0 

-9,7868 

352.3999 

3.0000 

48.0000 

0.0 

999.0000 

90.00 

7.1815 

35.00 

7.1815 

0.0 

0.0 

1 29 


0.0 

3.4052 

6.8104 

10.2157 

27.2418 

30.6470 

34.0522 

37.4574 

60.1989 

63.5641 

66.9693 

70.3746 

arr.4007 

90.8059 

94.21U 

97.6163 

-1.0783 

-0.8726 

-0.6739 

-0.5250 

-0.7616 

-0.9293 

-1.0423 

-1.1501 

-1.9552 

-2.0687 

-2. 1283 

-2.1737 

-1.8927 

-1.7281 

-1.5125 

-1.2571 

0.0 

16.0000 

30.2500 

45.0CCC 

100.0000 

100.0000 

too . 0000 

100.0000 

100.0000 

100.0000 

100.0000 

100.0000 

M.2900 

45.000C 

30.2500 

16.0000 

8«444446«4 

2 ♦44444 


94.4100 

6.5740 

-5.4960 

144.8940 

146.2500 

22.3940 

-5.4470 

92.5970 

4.0000 

10.0000 

1.0000 

0.0 

0.0 

0.0 

0 13 


0*0 

l.OOOO 

5.0000 

lO.OOOC 

60.0000 

70.0000 

80.0000 

90.0000 

0.0 

-0.0129 

-0. 1367 

-0.4893 

-3.1954 

-3.4439 

-3.6663 

-3.9C08 

0.0 

0.0529 

0.2700 

0.3078 

-1.6545 

-2.1469 

-2.6232 

-3.0692 


3 4444444444 


146.2500 

22.3940 

-5.4470 

92.597C 

176.8800 

31.7410 

-6.9000 

65.0C10 

3.0000 

10.0000 

1.0000 

O.C 

0.0 

0.0 

0 13 

0.0 

1.0000 

5 . 0000 

lo.ocno 

60.0000 

70.0000 

eo.oooo 

90.0000 

0.0 

0.0529 

0.2700 

0.3078 

-1.6545 

-2.1469 

-2.6232 

-3.0692 

0.0 

0.0138 

0.1706 

U.3492 

-0.3994 

-0.6585 

-0.9846 

-1.3231 




176.8800 

31.7410 

-6.9000 

65.0010 

190.5400 

36.8040 

-7.7160 

52,9«7C 

2.0000 

10.0000 

1.0000 

0.0 

0.0 

O.Q 

0 13 

0.0 

1.0000 

5.0000 

IC.OCOO 

60.0000 

TO.OOOC 

80.0000 

90.0000 

0.0 

0.0138 

0.1708 

0.3492 

-0*3554 

-0.6585 

-0.9846 

-1.3231 

0*0 

0.0151 

0.0774 

0.2623 

-0.2227 

-0.4511 

-0.7002 

-0.9606 


WIMG) MITH OAS NACELLES. (1/10 SCALEI. 

0.0 0.0 0.0 11 


0.0 1.000 
-5.4500 129.0000 


FUSELAGE (FUSIFORM PINED 
FUSELAGE (FUSIFORM PAHFLl 
NVORfRNCV, SPC.PDL 


-23.74 

7.1815 

-90.00 

7. 1815 


0 

1 


13.6209 

17.0261 

20.4313 

23.P365 

40.8627 

44,7679 

47.6731 

54.4835 

75.6044 

77.1850 

80.590? 

83.9955 

100.0000 




-0.4597 

-0.451? 

-0.4R81 

-0.5931 

-1.2789 

-1.4015 

-1.5233 

-1.7477 

-2.1793 

-2.1737 

-2.1254 

-2.0233 

-1.0499 




59.2900 

77.2500 

35.0000 

100.0000 

100.0000 

100.0000 

lOO.OOOC 

100.0000 

100.0000 

100.0000 

P^.GOOO 

7?.?5C0 


0.0 


WING PANEL NO. 1 
WING PANEL NH, 1 
NVORtRNCV , SPCfP^^L 
2 0 0 

20.00CC 30.0000 40,rf.OO 50.rooc 

LOO.OOOO 

-l.258« -1.9373 -2.4770 -2.884’ 

-4.1644 

0.1145 -0.2430 -0.6S25 -1.1609 

-3.4807 


WING PANEL 2 

WING PANEL NO. 2 

NVr»R*RNCVt SPCtPOL 
2 0 0 

20.0COO 30.C000 40.rCCO f>O.OCOO 

lOO.OCOO 

0.1145 -0.2430 -0.6&25 -1.1609 

-3. 4807 

0.4136 0.3277 C.1538 -C.C785 

-1.6677 


WING PANEL NO. 3 
WING PANEL NO. 
NVOR.RNCV.SPC.PPL 


0 

0 



20.0000 

100.0000 

30.00)0 

4C.C00C 

50.0000 

0.4138 

-1.6677 

0.3277 

0.1538 

-0.0785 

0.3095 

-1.2362 

0.2567 

0. 1373 

-0.0226 
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TABLE A-1 


Continued 



5 





190,5400 

36.8040 

-7.7160 

52.^^670 

WING PANEL 

NO. 4 


222.5800 

48.6820 

-8.9740 

31.1720 

WING PANEL 

NO. 4 


4.0000 

lO.OOCO 

1.0000 

0.0 

NVORtRNCV, 

SPCfPDL 


0.0 

0*0 

0 13 


2 0 

0 


0.0 

l.OOOC 

5.0000 

10.0000 

20.0000 20.0000 

40.0000 

50.0000 

60.0000 

70.00 00 

80.0000 

90.0000 

100.0000 



0.0 

0.0181 

0.0774 

0.2623 

0.3095 0.2567 

0. 1378 

-0.0226 

-0.2227 

-0.4511 

-0.7002 

-0.9606 

-1.2362 



0.0 

-0.0257 

-0. 1123 

-0.1 1*^5 

0.0513 0.1283 

C.1925 

0.2310 

0.2621 

0.29 19 

0,311? 

0.3?7-> 

0.3368 











222.5800 

48.6S20 -e.9740 


31.1720 


WING PANEL NO. 5 


248.0000 

64.5000 -lO.OOOC 


1^.4500 


WING PANEL NO. 5 


5.0000 

lO.OOOC 0.0 


c.o 


NVOR,RNCV,SPCtPDL 


0.0 

0.0 0 13 


2 

0 0 


0.0 

i.oono 5.0000 


10.C(*00 


20.0000 30.0000 40.0000 

50.0000 

60.0000 

70.0000 80,0000 


90.0000 


100.0000 


0.0 

-0.0257 -0.1123 


-0.1155 


C.0513 0.1283 0.1P93 

0.2310 

C.2631 

0.2951 0.311? 


0.327? 


0.3368 


0.0 

-0.0411 -0.1800 


-C.3137 


-0.4371 -0.4731 -0.4628 

-0.3651 

-0.2057 

-0.0051 0.2674 


0.5811 


0.9103 



7 






209.4000 

22.3940 -3.4000 


39.3000 


UPPER NACELLE PYLON 


209.4000 

22.3940 -5.4470 


39.3C00 


UPPER NACELLE PYLON 


1 .0000 

10.0000 0.0 


C.O 




0.0 

0.0 0 

0 


2 

0 0 









209.400 0 

22.39h0 3.3640 


39.3000 


UPPER NACELLE 


209.4000 

22.3990 3.3840 


39.3000 


UPPER NACELLE 


6 .0000 

lo.oocr 0.0 

999.0 coo 


NVnR,RNCV,SPC,PDL 


90.00 

3.392C 30.00 


3. 3920 


-30.00 3.3920 -90.00 

3.3920 

-150.00 

3.39‘>0 -’lO.OO 


3.3920 


-270.00 3.3920 


0.0 

0.0 1 

3 


2 

0 1 


0.0 

12.1522 100.0000 






0.0 

0.0 0,0 






100.0000 

lOO.OOOC 100,0000 







9 ***««««««* 






209. 1000 

22.3990 -5.9470 


39.4000 


LOWER NACELLE PYLON 


209.100 0 

22.3940 -7.2000 


39.4000 


LOWER NACELLE PYLON 


1.0000 

10.0000 0.0 


0.0 




o.u 

0.0 0 

0 


2 

0 0 



1C 






209.1000 

>2.3940 -7.2000 


39.4000 


LOWER NACELLE 


209.1000 

22.3940 -7.2000 


39.4000 


LOWER NACELLE 


6.0000 

lO.OOOC 0.0 


999.0000 


NVORfRNCVfSPCfPDL 


90.00 

3.2000 30.00 


3.2000 


-30.00 3.2000 -90.00 

3.2000 

-150.00 

3.2000 -210.00 


3.2000 


-270.00 3.2000 


0.0 

0*0 1 

3 


2 

0 1 


0.0 

11.9146 100.0000 






0.0 

0.0 0.0 






100.0000 

100.0000 100.0000 







ll 






248.0000 

64.5000 -10.0000 


19.45000 

FIN PANEL NO. & 


259.5398 

64.5000 -2.6000 


14.1100 


FIN PANEL NO. 1 


4.0000 

10.0000 0.0 


0.0 


NVORtR*<^*VtSPCfPDL 


0.0 

0.0 0 

0 


2 

0 0 



0 0 0 
<*•*«*•«* end of input deck 
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SUPERSONIC CRUISE VEHICLE <CL1607*3?/-r*4 wING) WITH PAT N 4CFLLES. <l/in <CALF^ 
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RECOMMENDATIONS FOR THE EFFICIENT USE OF THE VORLAX PROGRAM 


The definitions of the input and output data presented above should 
suffice for the running of any arbitrary configuration through the VORLAX 
program. But in order to achieve an efficient and accurate use of the capa- 
bilities provided by the generalized nonplanar vortex lattice method, the 
following recommendations should be adhered to: 

• When running a new configuration for the first time , run only one 
flow condition point, i.e. , NMACH=1 and NALPHA^l, to check out the 
input deck. This should be repeated for both subsonic and super- 
sonic Mach numbers if the analysis of the given configuration will 
extend through both flow regimes. 

• Use this relaxation solution (IS0LV=O) whenever possible. This is 
due to the fact that the relaxation procedure computational time 
varies somewhat proportionately to the second power of the total 
number of vortices making up the configuration, whereas the vector 
solution time varies as the third power. In addition, the vector 
solution involves an order of magnitude increase in the data trans- 
fer operations between the computer storage and core regions. 

• When major panels are in tandem, and lying approximately in the same 
plane, the spanwise distribution of vortex elements should be iden- 
tical to prevent the trailing vortex legs shed by the fore panel 
from running through the control points of the aft panel, in which 
case, pressure distributions may show spurious oscillations with 
possible consec^uent effects on aerodynamic coefficients. 

• In specifying the vortex lattice density of a given configuration 

as a whole, or of a particular major panel, by the quantities NV0R 
and RNCV, three different degrees may be considered: sparse, medium, 

and dense. To illustrate the case, the starboard panel of an aspect 
ratio 3 straight-tapered wing is sparsely latticed if NV0R=1O and 
RNCV=5, say; if NV0R=l6 and RNCV=10, then the panel may be said to 
have a medium density grid. A finely latticed, or dense, grid would 
be obtained for the values NV0R=3O and RNCV=20. In the above 
examples, the grid density was considered to be uniform in both the 
chordwise and spanwise directions. It is obvious that mixed-type 
lattices like dense-spanwise-sparse-chordwise are also possible. 

The type of lattice, when correctly specified, is a powerful tool 
for the program user. It allows him to achieve the lowest computa- 
tional cost for a given type of required data, as follows: 

(l) Stability and control type data (accuracy is essential only in 
the force and moment coefficients, good definition of the 
aerodynamic load distribution is not required, drag coefficient 
accuracy not too critical): sparse lattice. 
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(2) Load distribution and drag type data (good definition of the 
aerodynamic load distribution, both spanwise and chordvice, 
is required, accurate drag coefficient values are needed): 
medium density lattice. 

(3) Surface pressure distribution type data (good definition of the 

surface pressure distribution — control surface theory is 

essential, such as for wing design): lattice with maximum 

chordwise density and with medium spanwise density. 

• In addition to the lattice, or vortex grid, density, the type of 
vortex distribution (cosine and equal spacing) is highly significaiil 
in determining the accuracy of the results. It has been found that 
chordwise cosine spacing (LAX=0) is sui^erior in accuracy to tin- 
chordwise equal spacing (LAX=1). In terms of spanwise spacing, no 
significant differences have been obsei'ved between cosine (l,AY = 0) 
and equal spacing (LAY-l), the cosine spacing appearing to he 
slightly more accurate. For the same number of vortices, there 
is no difference in the computational cost between cosine and equal 
spacing lattices. 


• When the fusiform body representation is being used (PDL=360 and 
NPP=l), the number of sides of the polygon defining the body cros;; 
section should be kept as low as possible, e.g., the cross sectioi: 
of a body of resolution can be adequately represented by a hexagon. 
Also, when a very slender body with pointed nose is being consideron 
the nose (and afterbody if it is also pointed) should be arbitrarilv 
blunted in the input definition in order to minimize the numerical 
difficulties caused by the crowding of the vortical singularities in 
the body nose region. 

Figure A-7 illustrates a typical vortex lattice model of an advanced twin 
engine tactical fighter. This particular lattice model, with even chordwise 
and spanwise spacing (LAX=1, LAY=1), is considered sparse, and quite adequate 
for stability and control work, both longitudinal and lateral. In this model 
both the fuselage and nacelles are represented, or simulated, by flat plate 
elements, i.e., NPP=0. Obviously, this body simulation does not allow the 
computation of surface pressures, but it is adequate, and the most computa- 
tionally efficient, for stability and control work as well as for load dis- 
tribution and drag data. 




Figure A-7. - Vortex lattice panel representation for twin engine fighter 





APPENDIX B 


COMPLETE PROGRAM COMPILATION AND EXECUTION 

FOR 

A GENERALIZED VORTEX LATTICE METHOD 
FOR SUBSONIC AND SUPERSONIC FLOW APPLICATIONS 




HARDWARE AND SYSTEMS 


The VORLAX program has been run on several different IBM computer 
systems at Lockheed. All classified work conducted by the Advanced Develop- 
ment Projects Division, "Skunk Works", has been run on a 360/55 and no 
information can be presented about operation on this system. ^ All non- 
"Skunk Works" areas such as commercial and military engineering perform 
aerodynamic and loads analysis using VORLAX on the IBM 360/9I. Furthermore, 
a 370/168 has been used to assure that VORLAX can be efficiently run on this 
system if scientific applications are removed from the 360/9I system. 


The program has been used for a wide range of aircraft configurations 
and a good idea of the best operational procedures has been established. 

The program has been run at sizes ranging from d 5K words (260K bytes J o 
85K words dUOK bytes). An initial attempt was made to use a central memory 
region of 65K. This proved extremely inefficient because the number of 
accesses for l/O to peripheral storage soared with the imposed limitation 
of small buffers for each logical unit used. The program became I/O bound 
and larger buffers were called for. By going to large buffers, approximately 
3000 words are accessed with a single read or write. With the smaller 
buffers, the number of words per access might be 6OO - a factor of five 
smaller. When large cases are run on VORLAX, 5 to 10 million words may have 
to be accessed; it is readily apparent that the number of reads and writes 
can become overwhelming if the buffering is not adequate. 


Almost all cases run on VORLAX use less than 20 minutes of central pro- 
cessor time on the 360/91. Lockheed has demonstrated that a 25 minute case 
will run in about 15 minutes on a Control Data Cyber 175- It is, however, 
incumbent on the user to ensure that the system control language is properly 
adjusted to optimize operation of the program and the computer on which it 


is running. 

Since small central core size is usually an objective, it is worth noting 
that due consideration has been given to various methods of reducing VORW 
central memory requirements. All the obvious methods of core reduction have 
been found to have disadvantages. In particular, by studying the tree struc- 
ture of the program, it was determined that overlay would not significantly 
reduce size and it would introduce inefficiences . Consequently, in order to 
reduce core size requirements, the actual array sizes in the progr^ ave 
been minimized to allow reasonable configuration complexity without un- 
necessary waste. But more importantly, mathematical techniques were used 
in the problem solution that reduce the need to have large matrices 
central memory. One, two, or three elements are brought into central memory, 
processed, and shipped back to peripheral storage to minimize memory require- 
ments. This results in enormous savings in central memory. 


Since the Spring of 1975, the source coding has almost doubled in size 
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with a comparable increase in program flexibility and capability and yet 
central memory requirements have only increased about 20 percent. This was 
made possible only through continued attention to the optimization of the 
operation of the program on Lockheed's computer complex. It is of on-eoine 

Sfici^n^S possible means of Improving program operational 

ffr several novel ideas are being actively examined 

lor their possible future incorporation to the computer code. 
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PROGRAM 


UNIFIED VORTEX LATTICE METHOD 
FOR SUBSONIC AND SUPERSONIC FLOW 
(VORLAX) 

COMPLETE COMPILE 
AND EXECUTION 
IN CDC FORTRAN 
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C • • * 

C...CO» CLt CM» C.NJ 

c • • * 

C..-CX» CY 

C • • • 8L 
C • • • 

C...FN» FY 
IQ 

C*..NT 

C...RM 

C...SX 

c • • • xs 

C-t.YM 

C...YY 

c • • » 2C 
C...ZZ 

C • • *016 
C...COC* CNC 
C...CRMt CYM 
C • • • DCP 

C • • • 

C...EPS 

C...HAG 

C • » • 

C...ITS* JTS 
C...LAX*LAY 

C • • • 

C...NPP 

C...NXSf NYS* NZS 

c • • * 

C...POL 

c • • • 
c • • • 

C*..PSI 

c • • * 

C • • * PCS 

C*..RLH 

c • • • 

C..*SLE 

C...SPC 

c.*.tnl 
C • • * TNT 
C...VSP 
C...VSS 
C...VST 
C.. .XAF 

c • « * 

C • • • XTE 

C...ALFA 

C...ALOC 

C • • * 
c • • « 

C.*.8ETA 

c..*cbar 

C...CMTC 

c • • * 

C...CSUC 


force coefficients. 

PANFL X- AND Y*FORCE COEFFICIENTS. 

dihedral of chordwise wow of horseshoe vorti 

CES • 

PANEL NORMAL AND Y-FOWCE PER c 

angle of ATTACK AND MACH NUMBER INDICES. 

TOTAL NUMBER OF CHORDWiSE ROWS OF VORTICES. 

PANEL ROLLING MOMENT HEP UNIT Q. 

<%PANMISE LOCATION INDEX. 

abscissae OF FLO,! FIELD SURVEY CHOSS-PLaMES. 
PANEL YAVING MOMENT PER UNIT Q. 

y-coordinate of horseshoe vortex CENTROiOS. 
normal camber COOHOINATF in fraction of chord. 
Z-CoSrOINAtE OF horseshoe vortex CENTROIOS. 
MAXIMUM CHANGE PER 

rHOROWlSE CD AND CN I IMrS LOCAL CHOR )• 

PANEL ROLLING AND YAWING MOMENT COEFFICIENTS. 
PRESSURE coefficients (FTTHER LOAD OH 

accepta^e final relaxation cycle change, 
height above ground of moment reference 

panfT and strip FLOW EXPOSURE FLAGS. 

[MncrCHORDVlSE *NC SPANRISE DISTRIBUTION 

PANFI CHORDWISE NONPLANARITY FLAG. 

Sr S^FLOB FIELD SURVEY PLANES (ORThOGO 
NAL TO THE COORDINATE AXES). 

PANEL SPANWISE NONPLANARITY FLAG. ALSO PANEL 

pInIl dihedral if panel is flat in the SPAN 

YaS^ANGLE^^POSITIVE AHEN FLOB COMES from 
PORTSIOE) . 

CROSS-SECT ION RADIUS VECTOR. 

approximate value of eigenvalue radius of 
normalwash matrix. 

SLOPE AT LEADING EDGE. 

PANEL LEADING EDGE StCTlON FACTOR. 

TANGENT OF LEADING EDGE SWEEP. 

TANGENT OF TRAILING EDGE SWEEP. 

TANGENT OF SKEWED VORTEX LINE. 

SEMISPAN OF cHORnWiSE VORTEX ROW. 

actual horseshoe VORTEX rAMRFR 

CHORDWISE PERCENT COORDINATES AT WHICH rA BER 

TRAILING EDGF ABSCISSA OF VORTEX ROW. 
amai r OF ATTACK (IN RaOTANS). 

COMPONENT OF FREF-STREAM AND ONSET FLOWS NOR 
Mr?HHE LrFACE AT the CONTROL POINTS. ALSO 
AUXILIARY ARRAY. 

pranotl-glauept factor. 

CHORDWISE ?SnAL MOMFNT ABOUT OUARTEP CHORD 
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C • • aORAG 
C...HEAO 
C*..IOES 
C*.«IPAN 
C...ITER 
C...LIFT 
C...MACH 
C«..NPAN 
C«.«NVOR 
C...RNCV 
C • • • 

C...SLE1* SLE2 

C • • » 

C • « • SM AX 

C...SREF 

C...SURF 

C.*,VINF 

C*«.XBAR 

C**.XSUC 

C...YAWQ 

C...YNOT 

c • • » zbar 

C...2ETA 
C • • • 

C...ZLE1* ZLE2 

C • • • 

c • • • znot 

C«»«AINC1» A1NC2 
C*. .ALPHA 
C...COTOT 
C. • .CHORD 

C... 

C...CLTOT 
C...CMTOT 
C...CNTO'k 
C...CRTOT 
C...CYTOT 
C...DNDXlf 0 n 0X2 
C . * . 
c ♦ . * 
c . . « 

c *. .gamma 

c . . . 

C...ISOLV 

c . . . 

C...LESMP 
C...NMACH 
C... ONSET 
C . . . 

c . . . 

C...PSPAN 

c .. .rflag 

c . . . 

C...RNMAX 
C...ROLLQ 
C.. .SLOPE 
C.. .SYNTH 


PANEL INDUCED DRAG PER UNIT Q. 

PANEL description INFORMATION. 

DESIGN (SYNTHESIS) FLAG. 

HORSESHOE STRIP PANEL INDEX. 

ACTUAL NUMBER OF RELAXATION CYCLES. 

panel lift per unit Q. 

FREE-STREAM MACH NUMBER. 

NUMBER OF PANELS DEFINED IN THE DATA INPUT. 
NUMBER OF PANEL CHOROkISE STRIPS OF YOHTICES. 
CHORDWISE NUMBER OF HORSESHOE VORTICES FOR A 
GIVEN PANEL. 

SURFACE SLOPES AT LEADING EDGE OF PANEL SIDE 
EDGES. 

NT. 

CONFIGURATION REFERENCE AREA. 

PANEL SURFACE AREA (PLANFORM AREA). 
reference free STREAM VELOCITY. 

ABSCISSA OF moment REFERENCE CENTER. 
panel leading edge thrust PER UNIT 0. 
yam rate (0E6S/SEC). 

0UTTLINE origin OF FLOW FIELD SURVEY GRID. 
ORDINATE OF moment REFERENCE CENTER. 

INCIDENCE OF STREAMWISE VORTEX ROB (STRtP) 
CHOROLINE. 

leading edge offset of CAMBERLINE at THp 
panel SIDE EDGES. 

BATERLINE ORIGIN OF FLOW FIELD SURVEY GRID. 
CHOROLINE incidence AT SIDE EDGES OF PAnEL. 
angle of ATTACK (DECS). 

total induced drag coefficient. 

CHORD LENGTH MEASURED ALONG CENTRLINE Or 
STREAMBISE rob OF HORSESHOE VORTICES. 

total lift coefficient, 
total pitching moment coefficient, 
total yawing moment coefficient, 
total rolling moment coffficient. 
total side force coefficient. 

SURFACE SLOPES AT THE CONTROL POINTS MEASURED 
ALONG THE PANEL SIDE EDGES. IF DESIGN Is IN 
VOKEO* they are THE LOAD COEFFICIENTS A? THE 
LOAD points ALONG THE PaNEL SIDE EDGES. 
HORSESHOE VORTEX CIRCULATION . ALSO USED 
AS TEMPORARY DATA STORAGE. 

FLAG DETERMINING SYSTEM OF SOLUTION FOR THE 
BOUNDARY CONDITION EC3UATI0NS. 

PANEL LEADING EDGE SWEEP (OEGS) . 

NUMBER OF MACH NUMBERS PER CASE. 

ONSET FLOB VELOCITIES GeNERATEO BY THE 
ROTATION OF THE CONFIGURATION ABOUT THE 
MOMENT reference CENTER. 

PANEL SPAN. 

COEFFICIENT MULTIPLIER IN SYSTEM OF BOUnjOARY 
CONDITION EQUATIONS. 

number of VORTICES FOR A GIVEN CHORDBISp ROB. 
ROLL RATE (OEGS/SEO . 

SURFACE SLOPE AT THE LATTICE CONTROL POINTS. 
PANEL SYNTHESIS (DESIGN) FLAG. 
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nooon ooonooooonnnonooooo 


.TAPER 

.TITLE 

ucp AM 

.X apex I YAPFXt /APfK 
.CSTART 

.deltay* delta/ 

irtOATX* FLOATY 
. .INTRAC 
..INVERS 
. .lOUANT 
..ITRMAX 
..LATRAL 

..moment 

..NALPHA 
, .NPANAS 


..PHIMED 

IIpitchu 


^ANEL TAPER RATIO. 

CASE TITLE. 

configuration reference wing span. 

COORDINATES OF PANEL FlPST SIDE EDGE LEADING 
edge (PANEL APEX). 

CHORD length OF PANEL FIRST SIDE EOGF . 

RUFTLINE AND WATERLINE SPACING OF FLOW FIELD 
SURVEY GRID. 

VORTEX WAKE FLOATATION PARAMETERS. 

PANEL LATERAL NONPLANARiTY PARAMETER. 
design (SYNTHESIS) FLAG. 

PANEL SYMMETRY FLAG. 

MAXIMUM allowable NUMBER OF RELAXATION CYCLES. 
configuration OR FLIGHT CONDITION SYMMETRY 
FLAG. 

PANEL PITCHING MOMENT PER UNIT Q, 
number of ANGLES OF ATTACK PER CASE. 

to^al number of panels that have to he taken 

INTO account in THE ACTUAL COMPUTATION 
PROCESS. 

angular COOROINATE of cross-section RADtUS 
VECTOR. 

PITCH HATE (OEGS/SEC). 


|PER UNIT 01 means IPER UNIT FREE-STREAm 
dynamic PRESSURE I i 

c . • • 

c 
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74/74 0P1*1 


FTn 4.5^410a 


CONTROL *V«LX,AER0 

C « • • 

subroutine aero (EW* ITOTAL) 


•. .PURPOSE 


INPUT 


TO COMPUTE FORCt ANO MOMENT OAT A HY INTEGRATION nF 
pressure distribution data ANC By taking INTO ACCOUNT 

EDGE SIJCTION FORCFS. 


CALLING SE0U£NCE< 

EW s NOPMALWASH AT LEADING EDGE INFLUENCE COEFFiriFNT 
MATRIX (RETRIEVED RO'^ RY kOW FROM UNIT 9). 
COMMONS 

X, Y» 7. B2* DLt SXf YVt ? L% CCP, JTS* LAX* PSl* 

SPC* VSS* VST* XTE* ALFA* MACH* NPAN* NVOR* SREF. 
VINE* XRAR* YAWO* /bar* 7FTA* CHORD * GAMMA* LFSNp, 
PSPAN** RNMAX* ROLLO* SLOPE* TaPfR* CSTART* IOUAnT* 
LATRAL* NPANAS* pitcho. 


•OUTPUT 


calling sequences 

NONE. 


COMMON< 

CDf CL* CM* CN* CX* CY. Eg* 
CRM* CYM* CMTC* CSUC* DRAG* 
COTOT* CLTOT* CMTOT* CNTOT* 


FY* OM* YM* CnC* CNC, 
LIFT, SURF* XSUC* 
CRIOT* CYTOT* moment. 


. .SUBROUTINES 

..called none. 


..DISCUSSION SUBROUTINE AERO COMPUTES THE AERODYNAMIC FORCE AnD 
•• MOMENT COEFFICIENTS BY INTEGRAIik^G THE PRESSURE 

•• distribution and computing THE (.FADING EDGE SlJCTfON 

•• FORCES IN accordance WITH LAN/') PROCEDURE. THFRE ARE 

•• THREE C( ASSES OF COEFF IC I FNJS* As FOLLOWS < U) TOTAL 

.. CONFIGUPATION COEFFICIENTS. (2) PANEL COEFF I C T ENT S , 

•• and (3) 5TRIPWISE. OR CHORDWlSt* COEFFICIENTS. A| L 

•• total COEFFICIENTS ARE RFCERENCED TO WIND AXES. CLASS 

•• ANO < J) coefficients ARF REFfREnCEO either To RUDY 

•• OR To wind axes as RFUUIRFO by TmF corresponding COEFF, 

•• definition, aero subroutine is called by main FOp EVERY 

•• ANGLE OF ATTACK ANi) MACH NUMBEh; COMBINATION. 
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74/74 OPT=*l 


FTN 4.5^4lOtf 


10/04/76 lO.lS,!'^ 


CONTROL^VHLX .BOUNOY 

SUBROUTINE BOUNOY (I TOTAL) 

cIIIpurpose to calculate the onset floi« COMPONENI normal to 

ci;. THE BOUNDARY SURFACE AT THE VOkTfX LATTICE CON 

c!.. trol points. 


9/1 7/76 


C*. .PURPOSE 1 

C • • * 

C • • * 

C • • • 

c... INPUT ( 

C • • • 

c... ^ 

C » • • 

c • • • 

c * « • ^ 

c • • • 

c . ..output 

c . * . 

ciilsUBROUTINES 
C.. .CALLED 

C . . . 

cliioiSCUSSION 

C . . . 

c . . . 
c . • . 
c * • • 
c . . . 
c . . . 
c . . . 
c . . . 
c . . . 


CALLING 5EOUENCE < , ^ 

ITOTAL = TOTAL NUMBER OF HORSESHOE VORTICES. 

COMMON < 

ALFA* PSI* VINE* PITCHQ* ROLLC. YAWO. LAX* RNHAX, 
SX* CX* X* YY* ZZ* DL* CHORD* XHaR* ZBAR* SLOPE* 
RFLA6. 

COMMON < 

ALOC* onset. 


CALLED NONE. 

DISCUSSION THE ONSET FLO«rf COMPONENT NORMAL TO THE BOONOARY 

AT THE VORTEX LATTICE CONTROL POINTS IS CALCULATf 
BY projecting THE FREE-STREAM VELOCITY VECTOR 
ALONG The surface NORMAL AND TftKiNG INTO ACCOUNT 
A RIGID BODY ROTATION ABOUT T>t POINT (XBAR* 0* 
ZBAR). THE Onset flow normal CoMPMONEnT is ALUC. 
ONSET DENOTES THE RIGID BODY ROTATION INDUCED 
VELOCITY COMPONENT ALONG THE X-AXlS. BOTH ARRAY 
ALOC AND ONSET* ARE DIMENSIONLESS* I.E.* THEY aRF 
REFERENCED TO THE FREF-STREAM VELOCITY. 


747 

74H 

749 

750 

75 1 
75? 

753 

754 

755 

756 

757 

758 

759 

760 
7f>l 
76? 
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766 
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7f>H 
769 
77U 
771 
77? 

773 

774 

775 

776 

777 

778 

779 


109 



74/74 OPT=l EIN 4.S*A10a 

10/04/76 

10.15, 

CONTROL^VRLX 

.GAUSS «)-):J/76 

VORLAX 

856 

C • • • 


vorlax 

857 

subroutine gauss UTOTALt REAPER* fW* APT) 

VORLAX 

858 

C • • ■ 


VORLAX 

859 

C...PIIRPOSE 

TO SOLVE THE BOUNDARY CONniflOlsI rOUATIONS BY THE 

VORLAX 

860 

C • • • 

METHlt OF (CONTOLLEO SUCCESSIVE OVER-RELAXaT lONl , 

VORLAX 

861 

c • • • 


VORLAX 

862 

c.. . INPUT 

CALLING SEUUENCt< 

VORLAX 

86 3 

C • • • 

ITOTAL = TOTAL NUMBER OF HORStSHoE VORfJCES* 

VORLAX 

864 

c • • • 

REXPAR = RELAXATION PARAMETER. 

VORLAX 

865 

c • • • 

EW = ROW or NORMALWASh MATRIX. 

VORLAX 

866 

c » • • 

COMMON< 

VORLAX 

867 

c * * > 

CX. SX. LAX. ALOCf IDES. CHORD. RNMAX. INVERS. 

VORLAX 

868 

c • • - 

1 TRMAX . 

VORLAX 

869 

c • • • 


VORLAX 

870 

C.* .OUTPUT 

CALLING SEOUENCE< 

VORLAX 

871 

c • • * 

XRT = AUXILIARY VECTOR USEU IN C.S.O.R, SOLUTION, 

VORLAX 

872 

c • • • 

COMMON< 

VORLAX 

373 

c • • • 

BlG« DCP« EPS. HLM. SEE. ITER* GaMMA. SLOPE. 

VORLAX 

874 

c . • • 

GAMMA IS THE SOLUTION VECTOR OE BOUDARY CONDI TIOn 

VORLAX 

875 

c • • • 

EQUATIONS. I. E.. HORSESHOE VOpTEX CIRCULATION 

VORLAX 

876 

c • * . 

STRENGTHS. 

VORLAX 

877 

c • • • 

NOTES IF INVERS = 1 (OESIGN PROCESS) THEN GAMMA iS 

VORLAX 

878 

c * • • 

PART INPUT AND PART OUTPUT. 

VORLAX 

879 

c # . • 


VORLAX 

880 

C. . .SURKOUTlNFS 

vorlax 

881 

C. . .CALLED 

NONE . 

VORLAX 

88? 

C. * . 


VORLAX 

883 

C. . .DISCUS5T0N THIS SUBROUTINE SOLVES EOP THE CIRCULATION STRENGTH 

VORLAX 

884 

C . * • 

OF the horseshoe VORTICES THAT SATISFY THE B.C. OF 

VORLAX 

885 

c . * • 

NO MASS-FLUX ALONG THF NORMAL 10 THE SURFACE AT INF 

VORLAX 

886 

c . • • 

CONTROL POINTS. ThIS SOLUTION IS PERFORMED ITERA 

VORLAX 

887 

c . * . 

TIVELY ROW BY ROW BY USING THE C.S.O.R. METHOD. TF 

VORLAX 

888 

c * . . 

A GIVEN ROW IS PART OF A PANEL lo HE DESIGNED. 1,E,. 

VORLAX 

889 

c . • . 

IDES = 1. Then instead oe solving for gamma(=: xpr). 

VORLAX 

890 

C . • • 

The computation of the slope r i stribu j ion along 

VORLAX 

891 

c. . . 

that row is performed RY matrix mul t n l ICaT ton. 

VORLAX 

892 


SLOPE = Ew ^^GAMMA. 

VORLAX 

893 

c. . . 


VORLAX 

89 s 

C • . • 


VORLAX 

89S 
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74/7^ 0P1=1 


FfN 4.54410 a 


0/04/76 10.1S.13 


C0NTR0L*VRtX.6E0M 

SUBROUTINE GEOM (1 TOTAL! 

cIIIpurpose 

c « • • 


9/16/76 


TO COMPUTE The voptex lattice OEOMETRY ano surface 
SLOPES at the control POINTS. ALSO TO COMPUTE THf 
LOAD DISTRIBUTION AT THE CORRESPONDING LOAD POINTS 

IF design process is invoked <IDFS =1). 


c • • • 
c # . • 
c • * « 

c... INPUT 

c « • • 
c • • • 
c . • * 
c • • * 
c . . ■ 
c * . . 
c . • • 

c!!!output 

c . . . 
c . • * 
c * . . 
c » . « 
c . . . 
c • . * 
c . • • 
c . . . 

c.. .SUBROUTINES 
c «. .called none. 

c • • • 

c!I!oISCUSSION the VORTEX LATTICE GEOMETRY IS LAID OUT PANEL HY 


CALLING SEOUENCE< 

NONE. 

COMMONS 

DL* ITS. LAX, LAY, NPP, POL. PCS, VSS. NPAN, NVOw, 
RNCVt SLElf SLE2< ZLEl. Zl.E2» AlMCi* AINC2* DNOXi, 
0NUA2. GAMMA, LESWP. PSPAn* SYmTh. TAPER. XAPFX. 
YAPEX. /APEX. CSTART, INTrAC* lOuANT. LATRAL. 
PHIMEO. 

CALLING SEQUENCER 

ITOTAL = TOTAL NUMBER OF HORSESHOE VORTICES. 
COMMONS 

X. Y. Z. CX. DL< NT* SX. YY. ZZ. OCP* JTS. SLE« 
TNL* TNT* VSP. VSS* VST, XTE* ALOC* IPAN, SMAX, 
ZETA* CHORD* RNMAX* SLOPE, NPANAS. 

N0TE< DC AND VSS MAY BE EITHER INPUT OR OUTPUT 
DEPENDING ON CONFIGURATION CONDITIONS. 


C . . * 

c • • • 
c * . * 
c . . . 

C . . . 

c . . » 
c . . . 
c . . . 
c • . . 
c . * . 
c . . . 

C . * . 
C . • . 

c . . . 

C . . ■ 
C . . . 

c . . . 
c... 
c . . . 

C . . . 

c • . . 
c . * • 
c . . . 
c..* 

C a . . 
C • . • 


PANEL based on THE GEOMETRIC AND VORTEX OISTRIHU 
TION CHARACTERISTICS SPECIFIED FOR THE GIVEN PANfL 
IN the INPUT data UNPUT SUBROUTINE). EACH PANEL 
IS SUBDIVIDED INTO A NUMBpR OF X-AXIALWISE STRIPS 
(NVOR). EACH STRIP CONTAINING A GIVEN NUMBER tRN^4AX = 
RNCV) OF HORSESHOE VOrTICFS WEOSF BOUND TRAILING LFGS 
COINCIDE WITH THE X-AXIALWISE EDGES OF THE STRIP IF 
THE parameter NPP = 0. WHEN NPP .= 1 THERE IS NO 
longer a CONTINUOUS STRIP OF VORTTCES SINCE THEY APE 
not located IN THE SAME PLANE-* IN THIS CASE <NPP = 

D the STRIP BECOMES AN ARRAY OR ROW OF VORTICES 
LOCATED IN TANDEM BUT WHOSE SPANS ARE NOT NECES 
SAHILY EQUAL. EACH STRIP OR VOhTFX ROW IS IDENTI 
FIED BY AN INDEX (SX). EACH HORSpSHOE VORTEX IN A 
GIVEN STRIP OR ROW IS IDENTIFIED BY A SECOND INDrX 
(CX), THE VALUE CX = 1 DENOTING THE LEADING EDGE 
element, and CX * RNMAX = RNCV DENOTING THE LAST, 

OR trailing EDGE. HORSESHOE OF ThE ROW. THEREFORE 
EACH ANO EVERY HORSESHOE VORTEX yS UNIQUELY lOFN 
riEJE.O BY EITHFR AN OVERALL INoEX (WHICH RUNS FRoM 
1 TO ITOTAL) OR BY THF PAIR OF VALUES (CX, SX). 

THE SPATIAL LAY-OUT OF THF VORlEX LATTICE CDRRES 
PONDING TO A GIVEN PANEL OLPENUS ON THE VALUES Of 
TWO parameters^ POL AND NPP. IE POL .LE . OGO.O ImEN 
the transverse vortex segments of The same VALUE OF 

CX FORM A CONTINUOUS STRAIGHT LImE* BUT IF PUL 
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1079 

VORLAX 

1080 

VORLAX 

1081 

VORl AX 

1082 

VORLAX 

1083 

VORLAX 

1084 

VORLAX 

lOrtS 

VORLAX 

1086 

VORLAX 

1087 

VORLAX 

1088 

VORLAX 

1089 

VORLAX 
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VORLAX 

1091 

VORLAX 

1092 

VORLAX 

1093 

VORLAX 

1094 

VORLAX 

109S 

VORLAX 

1096 

VORLAX 

1097 

VORLAX 

1098 

VORLAX 

1099 

VORLAX 

1100 

VORLAX 

1101 

VORLAX 

1102 

VORLAX 

1103 

VORLAX 

1104 

VORLAX 

1105 

VORLAX 

1106 

VORLAX 

1107 

VORLAX 

1108 

VORLAX 

1109 

VORLAX 

1110 

VORLAX 

nil 

VORLAX 

1112 

VORLAX 

nil 

VORLAX 

U14 

VORLAX 

nis 

VORLAX 

1116 

VORLAX 

1117 

VORLAX 

1118 

VORLAX 

1119 

VORLAX 

1120 

VORLAX 

1121 

vorlaa 

112? 

VORLAX 

11?3 

VORLAX 

11 ?4 

VORLAX 

1 125 

VORLAX 

1 1?6 

VORLAX 

ll?7 

VORLAX 

1 l?8 

VORLAX 

1 129 

VORLAX 
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VORLAX 

1 111 

VORLAX 
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VORLAX 
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74//^ 0PT = 1 


FTN 


360,0 Then thf transvfpsf vomtex segments or samf rxt 
THOUGH still CONTINUOUS* FORM A POLYGONAL LINE 
WHEN PWuJECTED UN A PLANE NORMAL TO THE X-AXIS, IE 
NPP = 0 THEN ALL THE TRANSVERSE vORtEX SEGMENTS 
or A GIVEN ROW (SAME SXI LIE IN THE SAME PLANE* 

BUT IF NPP = I THEN ThF TpANSVERSE SEGMENTS OF A 
HOW ARF LAID ON THE ACTUAL BODY SURFACE. THE BOUNfO 
trailing legs or segments are always PARALLEL TO 
THE X-AxIS (UP TO THE TRAILING EoGE OF THE GIVEN 
STRIP OR ROW ) . 


10/04/76 

lO.lS 

VORLAX 

1134 

vorlax 

ins 

VORLAX 

1136 

VORLAX 

1137 

VORLAX 

1 I 3R 

VORLAX 

1139 

VORLAX 

1140 

VORLAX 

1141 

VORLAX 

1142 

VORLAX 

1143 

VORLAX 

1144 

VORLAX 

U4S 
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74/74 0P7=l 


riN 4.b^4lOA 


10/04/76 10. IS. 13 


C0NTR0L*V«LX.MAP 

subroutine map Hit* EWXf EWY* lTOT/yL> 
c!!!purpose to compute The flow field abolf the configuration. 


c... INPUT CALLl 

FW * 

c • • • 

C.,. EWX = 

c • • • 

C... EWY > 

c » • • 

C... ITOT^ 

c... COMM( 

C.., IH* 1 

C... ZN0T« 

cil.OUTPUT CALL 

c... none 

C,.. C0MM( 

€••• NONE 

C... DIRE( 

c . • • X s ♦ 

c • » « 

c,.. vx» ' 

c • • • 
c • • • 
c » • ■ 

c.,. EPSLI 

C... SIGH 

C,.. CP » 

c.., RH « 

c... PPTO 

C... PIE 

c • • • 

c. ••SUBROUTINES 
€•• •CALLED NONE 


CALLING SEQUENCES 

EW * UPWASH INFLUENCE COEFFICIENT MATRIX (REIRIEwED 
ROW rY row from unit 3) . 

EWX = AXIALWASH influence COEFFTCTENT matrix (RETRIEVED 
ROM BY ROW FROM UNIT 4>^ 

EWY * STDEWASH INFLUENCE COEFFICIENT MATRIX (RETRIEVED 
ROW BY ROW FROM UNIT ?)• 

ITOTAL = total NUMBER OF HORSESHOE VORTICES^ 

COMMONS 

IH* IQ, NT* XS* NXS, NYS* NZS* PSl* ALFA* MACH* VNOT, 
ZNOT* alpha* GAMMA* RNMAX* DELTAy* OELTAZ. 

CALLING SEQUENCES 
NONE. 

COMMONS 

none. 

DIRECT PRINTS 

XS* YKl, ZKl = FIELD GRID POIM COORDINATES (HODv AXIS 
SYSTEM) . 

VX* VF* WF = TOTAL DIMENSIONLESS (REFERENCED TO VEI.O 

city AT iJPST'^EAM INFINITY) VELOCITY COMPO- 
NENTS along the x-y-z axes respectively 
(BODY AXIS SYSTEM) . 

epslon = uPwASH Flow angle in degrees. 

SIGMA = SIOEWASH FLOW ANGLE IN DEGREES. 

CP » LOCAL PRESSURE COEFFICIEM. 

RM » LOCAL MACH NUMBER. 

PPT07 = (LOCAL STATIC PRESSURE >/( TOT AL PRESSURE). 

PIF * (LOCAL STATIC PRESSURE )/< FREE STREAM STATIC PRES 
SURF) . 


..DISCUSSION FLOW FIELD QUANTITIES ARE COMPUlpO AT THE NODAL JOINTS 
OF A 3-O-GRlO DEFINED AROUND THE CONFIGURATION Hv A SET 
OF orthogonal planes ( X = CONST.* Y = CONST.* AnD 
Z « const. PLANES). ThE VFLOCIUfS ARE CALCULATEn nY 

the use OF Influence coefficient matrices based dn the 
vortex lattice represent ATION uF THE configuration. 
THESE matrices ARE COMPUTED IN SUBROUTINE SURVEY AND 
STORED IN UNITS 3. 4* AND 7 (OnE MATRIX AND ONE UNIT 
.. per velocity COMPONETl. ThE PRESSURE RATIOS AND RELATED 

FLOW QUANTITIES ARE COMPUTED THROUGH THE USE OF tSFN 
.. TROPIC FLOW relationships. 
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1996 

1997 

1998 

1999 

2000 
2001 
2002 
20 0 3 
2004 

2 ons 
2006 
20 0 7 
2008 
20 09 
2010 
2011 
2012 
201 3 
20 14 
20 IS 
20 16 

2017 

2018 

2019 

2020 
2021 
2022 
2023 


113 



74/7A OPT*l FTN 4.5^410 a 

10/04/76 

10.15.13 

CONTROL*VRLX 

•MATRX 9/16/76 

VORLAX 

2147 

C • • • 


VORLAX 

2148 

SUBROUTINE MATRK (EW* EUt ITOTAL) 

VORLAX 

2U9 

C • • • * 


VORLAX 

2150 

c... PURPOSE 

TO CiENEHATt THREE AERODYNAMIC INFLUENCE COEFF. MaTRI 

VORLAX 

2151 

c • • « 

CESS (1> THE NORMALWASH AT THE CONTROL POINTS 

VORLAX 

2152 

c • • • 

/ UNIT (2) THE AXIAL wash AT THE CONTROL POINTS 

VORLAX 

2153 

c • • • 

(EU / UNIT 2)4- AND (3) THF NOfiMAf WASH AT THE LEAo. EDGE 

VORLAX 

2154 

c ■ « • 

(EW / Unit 9), THESE MATRICES kEpRESENT THE INDUCED 

VORLAX 

2155 

c • • • 

velocity field due to THE HORStSHOE VORTICES OF THE 

VORLAX 

2156 

c » • • 

LATTICE. 

VORLAX 

2157 

c... 


VORLAX 

2158 

C.. .INPUT 

CALLING SEQUENCES 

VORLAX 

2159 

c . • • 

ITOTAL = TOTAL NUMBER OF HORSESHOE VORTICES. 

VORLAX 

2160 

c • • . 

COMMONS 

VORLAX 

2161 

c • . . 

Xf H2f rX* OL* NT. SX. YY. Z/. HAG* LAX. NPP* PSf. TNT. 

VORLAX 

2162 

c • . . 

VSP. VST. ATE. ALFA* IPAN. XHAK. ZbAR. CHORD. RF| AG* 

VORLAX 

216 3 

c • • « 

RNMAX. slope. FLOATX. FLOATY. INvERS. LATRAL. 

VORLAX 

2164 

c..« 


VORLAX 

2165 

C.. .OUTPUT 

CALLING SEQUENCES 

VORLAX 

2166 

c . . . 

EW = COigTROL POINT NORMALRASH MATRIX (STORED ROW BY 

VORLAX 

2167 

c . . . 

ROW IN UNIT 1 > . 

VORLAX 

2168 

c • . . 

EW * LEADING EDGE NORMALWASH MATRIX (STORED ROW 

VORLAX 

2169 

c . . . 

ROW IN UNIT R) . 

VORLAX 

2170 

c . . . 

EU •* AXfALWASH MATRIX (STOKED ROW BY ROW IN UNIT 2). 

VORLAX 

2171 

c . . . 

COrtMON< 

VORLAX 

2172 

c . • • 

NONE. 

VORLAX 

2173 

c . . . 


VORLAX 

2174 

C.. .SUBROUTINES 

VORLAX 

2175 

C.. .CALLED 

WASH. UXVEL. 

VORLAX 

2176 

C . . . 


VORLAX 

2177 

C. ..DISCUSSION THE ELEMENTS OF THt INFLUfncE COEFFICIENT MATRICfS 

VORLAX 

2178 

C . . . 

ARE generated BY cOMPirTlNG THE CORRESPOND TNG VELOCITY 

VORLAX 

2179 

C . . . 

INDUCED AT The (Kl.Jl) CONTROL POINT BY THE IK.J) 

VORLAX 

2180 

C . . . 

HORSESHOE VORTEX OF UNIT STRENGTH. IF K = K1 AND 

VORLAX 

2181 

c . . . 

j a Ji (SElf-influencf j Then the principal part of 

VORLAX 

2182 

c . . . 

The oownwash Integral is aooec to the computatiom or 

VORLAX 

2183 

c... 

the corresponding EW coefficient, also if the COmTROL 

VORLAX 

2184 

c . . . 

POINT IS within a given NEAR FjETt D RADIUS OF ThE 

VORLAX 

21HS 

c . . . 

INDUCING HORSEShOF VORTEX. THE AXIALWASH C(Jrg T R I BoT T ON 

VORLAX 

2186 

c . « . 

IS COMPUTED BY IN TERD I G I T A I EO VORTEX SPLITTING. 

VORLAX 

2187 

c . . . 


VORLAX 

2188 

c . . . 


VORLAX 

2189 



74/74 0PT=1 


FTN 4,5*410a 


10/04/76 10.15.13 


CONTROL^VRLX. input 

SUBROUTINE INPUT 

c • • • 

c!!. PURPOSE TO read 

C... USE IN 


C • • • 

c • • • 

c... INPUT 

c . • • 
c • • • 
c • • • 
c . . . 

C.. .OUTPUT 

c . * • 
c * * * 
c * * * 
c » * * 
c . * . 
c • • * 
c . * « 
c . . • 
c . • • 
c • • • 


TO read in INPUT data AND PREPaRF SUCH DATA FOR 
USE IN THE GENERATION OF VORTEX LATTICE GEOHETRV 
TO BE DONE IN SUBROUTINE GEOM. 


CALLING SEOUENCE< 

NONE. 

COMMONS 

LAX. 

CALLING SEQUENCES 
NONE. 

COMMONS 

DL* XS. ITS. NPP* NX5. NYS* NZS. POL* PSI* RCS* 

SPC* vss* C0AR. Head. mach. npan. nvor. rncv. slfi. 
SLE2. SREF. VINF* XBAR. YAWO. YNnT. ZrtAR. ZLEi. 
ZLE2. ZNOT. AINCI. AINC2* ALPHA* ONDXl. DNDX2. 
gamma. I ESWP. NMACH. P5PAN. RCLLO. SYNTH. TAPER* 
WSPAN. XAPEX. YAPEX. ZAPEX. CSfAPT* UELTAY. DELTaZ. 
INTRAC. INVERS. lOUANT. LATRAL. nALPhA. PHlMEO* 
PITCHQ. 


cliisUBrtOUTINES 
C... CALLED NONE. 

c!I!oISCUSSION a master frame of reference is ASSUMED IN DEFINING 
A CONFIGURATION. THIS FRAME OF WEEERENCE IS AN 
r... orthogonal CARTESIAN COORDINATE SYSTEM, THE X-Z 

r.*. BEING THE CENTERLTNE PLANE H M THE X-AXIS POINliNG 

doi^nstrfam. and The z-axis directeu uphiaro*^ the 

Y-AXIS POINTS TO STARBOARD. THE ORIGIN OF THE SYs 

C,,, tem can be any convevient point in the x-z plane. 

r... The configuration can be mde up of symmetrical 

r* (about the x-z PLANE) AND/OR ASYMMETRICAL COMPO- 

C*.. NENTS, and in defining THF SYMMETRICAL COMPONENTS 

only THF starboard ELEMENTS ARt SPECIFIED. THE 
C... CONFIGURATION TO RE INPUT IS DIVIDED INTO A SET 

OF MAJOR PANELS* UP TO 20 OF IhESF PANELS CAN RE 
r... input, SYMMETRICAL COMPONENTS (LFFT ♦ RIGHT) HEImG 

r... COUNTED ONLY ONCE. FOR INSTANCE* A WING OF ZERO 

r... THICKNESS AND WITH STRAIGHT LEADING AND TRAILING 

r... EDGES, AND WITH LINEAR LOFTING BFTWEEN ROOT AND 

r TIP* CONSTITUTES A MAJOR PANEL. COMPLEX PLANFORH^;, 

^*** ANO NON-LINEAR CHANGES IN TWIST AND AIRFOIL SFCTiOMS 

are DESCRIBED BY DEFINING MORE THAN ONE PANEL FOk- A 
-*** given wing. SUBROUTINE INPUT PREPARES ThE UATA SuF 

. CIFIED FOR EACH MAJOR PANEL SO ThAT THEY CAN LATpr 

be used IN SUBROUTINE GEO 4 TO oENERATE THE PWOPt^J 
rl* VORTEX I ATTICE FOR EACH PANEL. 

r AN AIRFOIL WITH THICKNESS CAN HE REPRESENTED BY a 

P*** double vortex Sheet, i. e.. by defining two majou> pa 

P* nels arranged in a j«RIPLAnE^*. or ^sanowich^, fashion* 

r .. one panel representing THE UPPER SURFACE OF THE 

C FOIL, AND the other PANEL REPRESfNTING The LOwFR 
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C.* 

C • • 

c • • 



c.. 
c*. 
c.. 
c • • 
c.. 

C • 

c.. 
c.. 
c « • 
c • • 
c.. 
c • • 
c ■ • 


FTN 4.S+410A 


SURFACF OF THE SAME AIRFOIL. 

fusiform bodies are modelled Bt defining an ADXIi iary 
BOUYIOENTICAL in CKOSS-SEC f lOKAL SHAPE AND LONGlruoI 
NAL area distribution TO THE ACTUAL BOOT BUT ^IIhOUI 
CAMBER, the AUXILIARY BODY CROSS-SECTION IS APPROXI 
MATED RY A P0LYG0(4 ii^HDSE SIDES DETERMINE THE IRAmS 
VERSE LEGS OF THE HORSESHOE VORTjCES. THE VERTICES 
OF THE POLYGON AND THE AUXILIARY BODY AXIS DEFINf 
A SET OF RADIAL PLANES IN WHICH THI' HOUND IWAlLUtG 
LEGS OF THE HORSESHOE VORTICES LIE PARALLEL To ThE 
AXIS, AS THE CROSS-SECT IOn CHAnGfS SHAPE ALONG IhE 
AXlSf THE polygon CHANGES ACCORDINGLY HUT WITH ThE 

constraint that The Polygonal vertices must always 

LIE IN the SAME SET OF RAoTaL pLaNES. THE BOOr CaMmER 
IS spectfieo inuependfntly. 
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10. IS 

VORLAX 

1473 

vorlax 

1474 

VORLAX 

147S 

VORLAX 

1476 

VORLAX 

1477 

VORLAX 

1478 

VORLAX 

1479 

VORLAX 

1480 

VORLAX 

14Hl 

VORLAX 

14H? 

VORLAX 

148 3 

VORLAX 
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148S 
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\0/0**/lt> 10«1S»13 


C0NTR0L*VRLX .PRESS 

subroutine press UTOTALf EU) 


4/15/76 


• .PURPOSE TO COMPUTE PRESSURE LOAD COEFFICIENTS (CPLOl^ER 

CPUPPER)* OR SURFACE PRESSURE COEFFICIENTS* FROM THE 

VALUES OF the Induced velocities and circulation 
.. strengths. 


.•INPUT 


, .OUTPUT 


CALLING SEOUENC£< 

ITOTAL = TOTAL NUMBER OF HORSESHOE VORTICES. 

EU a AXIALVASH influence COEFFICIENT MATRIX (RETRIEVED 
ROW RY ROW FROM UNIT 21. 

COMMONS 

B2* CX* DL* SX* JTS* PSI* TNL* TnT* ALFA* YAWU* rHORD. 
GAMMA* ONSET* RNMAX* SLOPE* WSPAN. 

CALLING sequencer 
NONE. 

COMMONS 

OCP. 


.. .subroutines 

...CALLED NONE. 


!!!oISCUSSION this SUhROUTINE COMPUTES THE PRESSURE COEFFICIENT ARRAY 


OCP. EACH ELEMENT OF DCP CORRESPONDS TO A GIVEN 
horseshoe VORTEX OF THE LATTICE* AND IT IS ASSUMpO TO 
ACT AT THE HORSESHOE CENTROID* I.E.* THE MI= 

THE transverse* OR SKEWED. LEG. AN ELEMENT OF THF DCP 
ARRAY IS either A LOAD COEFFICIENT (IF THE SURFACE 
IS ASSUMED WETTED ON BOTH FACES. I.E.* JTS = ITS = D* 
OR A SURFACE PRESSURE COEFFICIENT (IF THE SURFACE IS 
ASSUMMED WETTED ON ONE SIDE ONLY. I.E.* JTS = iTs = 0). 
IN THE COMPUTATION OF LOAD COEFFICIENTS* THE ClRcULA 
TION strengths* The free stream VELOCITY components* 
and the onset flow due to any angular rotation Of the 

CONFIGURATION ARE TAKEN INTO ACCOUNT^ THE EFFECT OF 
THE INDUCED AXIALWASH IS IGNORED (SECOND ORDER ERROR). 
IN the COMPUTATION OF SURFACE PRESSURE COEFFICIENTS* 

The axialwash is also taken into account, if npp = i* 
the #axiaukasha is no longer a True axialwash <x-axis 

VELOCITY COMPONENT)* RATHER IT REPRESENTS THE VEiOClTY 
component tangential to the surface BUT WITH THE SIDE 
WASH left out. surface PRESSURE COEFFICIENTS ARE LIMI 
TED BY THE ISENTROPIC VALUES CORRESPONDING TO SiflGNA 
TION and 70 percent of VACUUM FOR THE GIVEN FREE-STREAM 
MACH number. 
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7A/7<fr 0PT»1 FTN ^.5+410 a 

lO/OA/76 

10.15 

CONT«OL*VRtX 

•PRINT 9/17/76 

VORLAX 

2768 

C • • • 


VORLAX 

2769 

SUBROUTINE PRINT (ITOTAL* EW* FWXt EWY) 

VORLAX 

2770 

C • • • 


VORLAX 

2771 

C... PURPOSE 

TO PRINT program t)AlA, 

VORLAX 

2772 

c 


VORLAX 

2773 

C*.. INPUT 

CALLING sequencer 

VORLAX 

277A 

C • • • 

ITOTAL = TOTAL NUMBER OF HORSESHOE VORTICES. 

VORLAX 

2775 

c • • • 

EWt EWX. E<«*Y = STORAGE ARRAYS SET ASIDE FOR CALLjNG 

VORLAX 

2776 

c • • • 

SUBROUTINE MAP. 

VORLAX 

2777 

c • • « 

COMMONS 

VORLAX 

2778 

c • • • 

CD# CL* CM* CM* CX* CY* DL * IH* TO* YY* ZC* ZZ* 

VORLAX 

2779 

c • • « 

BIG* COC* CNC* CRM* CYM, OCP* EPS* HAG* JTS* LAX. 

VORLAX 

2780 

c • • • 

NXS* POL* P5I* RLM* SPC* CBA«* CmTC* CSUC* HEAD* 

VORLAX 

2781 

c • • « 

IRAN* ITER* LIFT* MACH* NOAN* NVnR* RNCV* SREF* 

VORLAX 

2782 

c • • « 

SURF* VINF* X0AR* XSUC* YaWO* ZBAR* AlNCl* AINC?* 

VORLAX 

2783 

c • • • 

ALPHA* CDTOT* CHORD* CLTOT* CMfOT* CRTOT* CYTOT* 

VORLAX 

2784 

c • • • 

ONOXl* DNDX2* GAMMA* ISOLv* LE:>WP* PSPAN* RNMAX* 

VORLAX 

2785 

c • • • 

ROLLO* SLOPE* SYNTH* TAPER* TITLF* WSPAN* XAPEX* 

VORLAX 

2786 

c • • • 

YAPEX* ZAPEX* CSTaRT* FLOaTX* FLOATY* INVERS* 

VORLAX 

2787 

c • • • 

IQUANT* IIRMAX* LATRAL* MOMENT* PITCHO. 

VORLAX 

2788 

c • • • 


VORLAX 

2789 

C*.. OUTPUT 

CALLING sequences 

VORLAX 

2790 

c • • • 

NONE* 

VORLAX 

2791 

c • • • 

COMMONS 

VORLAX 

2792 

o« • • 

NONE* 

VORLAX 

2793 

c • • • 


VORLAX 

2794 

C*.. SUBROUTINES 

VORLAX 

2795 

C... CALLED 

ZNORH* MAP. 

VORLAX 

2796 

C • • • 


VORLAX 

2797 

C. . •DISCUSSION THE DATA PRINTED OUT HY ThIS SUBpOUUNE ARE ARRANGED 

VORLAX 

2798 

C • « » 

IN three groups* as followss 

VORLAX 

2799 

c • • • 

(1) INPUT Data* these are data which are taken From 

VORLAX 

2800 

c • • * 

the input DECK ANO HAVE BEEN CONVERTED IN A FORMAT 

VORLAX 

2001 

c • • • 

SUITABLE FOR USE BY THE PR(JGpAM* 

VORLAX 

2802 

c • • • 

(2) PANEL AND TOTAL CONFIGURATION FORCE ANO MOMENT 

VORLAX 

2803 

c • • • 

COEFFICIENT DATA. 

VORLAX 

2804 

c ■ • • 

(3) DATA WHICH ARE RELATED EIThEr TO CHOROWISE STRIPS 

VORLAX 

2805 

c • • • 

OR TO INDIVIDUAL HORSESHOE VORTICES. 

VORLAX 

2806 

c • • • 

THESE THREE DATA GROUPS ARE PRINTED SEQUENTIALLY. 

VORLAX 

2807 

c • • • 

IF A FLOW FIELD SURVEY HAS BEEN REQUESTED* THE CoRRES 

VORLAX 

2808 

c • • • 

PONDING DATA ARE PRINTED OUT THROUGH SUBROUTINE MAP. 

VORLAX 

2809 

c • • • 


VORLAX 

2810 


U8 



7A//^ 


OPT - I 


f T 4 • 5 < 4 1 0 A 


10/04/76 10. IS. 13 


CONTROL ‘^VKLX.SUWVeV 

C. . . 

SiiHkOUfINL SUkVfCY (LW. f LV /^ IIUTAD 


^/17/76 


r. . 
c. . 

c. , 

r, 
c. 
c. 
c. 
c. 
c . 
r, . 
c. 
c. 
c. 
c. 
c. 
c. 
c. 
c. 
c. 
c. 
c. 
c. 
c. 
r . 
c. 
c. 
c. 
c. 
c . 
r. 
c • 
c. 
c. 
c • 
c. 
c. 
c. 
c. 
c • 
c. 
c • 
c. 
c. 


PUHPOSe TO OtNEkArt THkEF: AtttOOYNAMiC INFLUENCE COfcFFlClFNT 
MATkICES< (U the UPWASH aT Tht FLOW FIELD SURVEY 
POINTS (FW / unit 3)^^ <2) THE a^tALWASH AT THE F| OW 
FIELD SURVEY POINTS (FWK / UN 1 I 4 ) '*' AND (1> THE 
SIDEWASH aT the Fl OW FIELD SURVEY POINTS (EWY / UNU 7). 
THESE MATRICES REPRESENT THE INDUCED VELOCITY FlFLD 
DUE To THE HORSESHOE VORffCES OF THE LATTICE. THiS FLOW 
FIELD IS MtASURED AT THE NODAL POINTS OF A SPECIFIED 
3-U GRID. 

input CALLING SEUUENCt< 

ITOTAL -■ TOTAL NUMBE^-^ OF HORSE SHOE VOwTICES. 

COMMON< 

X, B^• CX* DL» NT, SX, XS. YY* Z7« HAG* LAX, NXS, 

, NYS* NZS, PSI, TMT, VSP« VST* XTp* ALFA» YNOT , XhAR* 

. /BAR* ZnOT, chord* RNMAX* i)ELTay* deliaz* FLOATX* 

floaty* IATRAL* 
iourpur calling sfoofncex 

EW = UPwASH MATRIX (STORED ROW HY ROW IN UNIT 3) 

’ tWX = AXIALWASH matrix (STORED ROW BY ROW IN UNIT 4) 

[ EWY = SIDEWASH MATRIX (STORED ROW HY ROW IN UNIT 7) 

, COMMON < 

NONE* 


*SUHROUTlNt S 

.CALLED WASH* UXVEL. 

[discussion the EI.EmENTS of The ImFlUknCE LOKFFICIENT MATRICtS are 
generated by computing IMF VELOCITY INDUCED AT A F(,OW 
FIELD SURVEY POINT HY THE (K* J> HORSESHOE VORTEx OF 
UNIT strength, if THE FLOW FIELD SURVEY POINT 1= 

WITHIN A GIVEN NEAR FIELD RADIUS OF T HF. INDUCING 
horseshoe vortex THEN THE AXIAlWaSH CONTRIBUTION IS 
COMPUTED HY 1 NTERD I G I T A TFn VOMEx SPLITTING. THE Ft OW 
FIELD SURVEY POINTS ARE IHE NOOAj POINTS OK A 3*0 GRID 
DEFINED BY A SET OF THREE ORlhi)(iONAL PLANES, THESE 
PLANES ARE SPECIFIED HY THE INPUT VALUES OF XS* YNOT ♦ 
OELTAY* ZNOT* AND DEt lAZ. 
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10/04/76 

10.15.13 

C0NT«0L<»V«LX 

•UXVEL 9/14/76 

VORLAX 

3355 

C * • • 


VORLAX 

3156 

subroutine uxvel z* S* T» R?* TOlZ* U) 

VORLAX 

3357 

c • • • 


VORLAX 

3 35B 

c.. .PURPOSE 

TO COMPUTE The axialwash tnouceo by a skewed 

VORLAX 

3359 

c . • • 

RECTILINEAR VORTEX SEGMENT OF UNIT CIRCULATION. 

VORLAX 

3360 

c • • » 


VORLAX 

3361 

C - . . INPUT 

CALLING SEOUENCES 

VORLAX 

3362 

c « * • 

X. Y* 7 = ORTHOGONAL CARTESIAN COORDINATES OF 

VORLAX 

3363 

c • • • 

receiving point measured in a referenck frame 

VORLAX 

1364 

c • • • 

CENTERED AT- THE MIDPOINT OF VORTEX SFG^EMT^- 

VORLAX 

3365 

c. .. 

THE X-AXIS IS PARALLEL TO THE X-AXIS Op THE 

VORLAX 

3 366 

c • . • 

MASTER (CONFIGURATION) COORDINATE SYSTpM. THE 

VORLAX 

3367 

c. .. 

Y-axis is normal to thf x-axis hut lies in 

VORLAX 

3368 

c • • • 

The Plane determined by the x-axis and the 

VORLAX 

3369 

c • • • 

VORTEX SEGMENT. AND THF Z-AXIS IS NORMAL TO 

VORLAX 

3370 

c • • • 

SUCH PLANE. 

VORLAX 

3371 

c . • • 

S » SEMlSPAN OF VORTEX SEGMENT. 

VORLAX 

337? 

c • • • 

T = tangent of sweep angle of vortex segment. 

VORLAX 

3373 

C • • • 

02 a COMPRESSIBILITY FACTOR (» MaCH **2 - l.O), 

VORLAX 

3 374 

c* . . 

TOLZ = ^^UMtRICAL TOLERANCE CONSTANT. 

VORLAX 

3375 

c • • • 

COMMONS 

VORLAX 

3176 

c • • * 

none. 

VORLAX 

3)77 

c • • • 


VORLAX 

3378 

C... OUTPUT 

calling seuuences 

VORLAX 

3379 

c • • • 

U * X-AXiS velocity component (AXTALWASH) induced by 

VORLAX 

33B0 

c • • • 

SKEWED VORTEX SEHMENT OF UNIT INTENSITY. 

VORLAX 

3 381 

c. . • 

COMMONS 

VORLAX 

3382 

c • • • 

NONE. 

VORLAX 

3383 

c • • • 


VORLAX 

3384 

C- . .SUBROUTINES 

VORLAX 

3385 

C. . .CALLED 

NONE • 

VORLAX 

3386 

c • * • 


VORLAX 

3387 

C.. .DISCUSSION THIS SUBROUTINE IS CALLED WHEN ThE AXIALWASH IS 

VORLAX 

3188 

C « • • 

COMPUTED IN accordance WITH VORTfX SPLITTING SCHrMF. 

VORLAX 

3389 

c . . . 

ONLY THF. AXIALWASh IS COMPUTED. ANO ONLY THE TRAvS 

VORLAX 

3 390 

c • • . 

VERSE SEGMENT OF THE HORSESHOE Is TAKEN INTO ACCOUNT. 

VORLAX 

3391 

c . . . 

THE SAME COMMENTS PRESENTED IN SUBROUTINE WASH 

VORLAX 

3 392 

c • . • 

REGARDING THE NUMERICAL SINGULARITY IN THE VICINTTy 

VORLAX 

3393 

c . . . 

OF THE CHARACTERISTIC SURFACES (MACH CONES) ARE 

VORLAX 

3394 

c . . • 

ALSU APPLICABLE HERE. SUHR0UTIr4E UXVEL IS A OIRErl 

VORLAX 

3395 

c . . . 

COPY or parts of subroutine WASH. 

VORLAX 

1396 

c . . . 


VORLAX 

3397 

c . . . 


VORLAX 

3398 


120 
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CONTROL*VRLX. VECTOR 9/14/76 

c • • » 

subroutine vector (ITOTALt NAl. AV. EW# VORl* VORK* VOkD 
C 
C 

c. . .PURPOSE 

c . . . 
c.. . 
c . « • 

c.. . INPUT 

c . * . 
c . . . 
c . . . 
c . . . 
c . . . 
c . . . 

c. . .OUTPUT 

c . > . 
c . • • 
c . . . 
c . . . 
c . . . 
c . . . 

c. - .subroutines 

C... CALLED NONE. 

C. . .DISCUSSION the BOUNDARY CONDITION EOtiATlOnS ARE SOLVED DIRECTLY 

C... BY A VECTOR OR FMOtiONAL I2AT ION PROCEDURE (PURCELL/S 

C... VECTOR METHOD). SETS OF LINEARLY INDEPENDENT VECTORS 

C... ARE constructed WHICH ARE SUCCtSlVELY ORTHOGONAL TO 

C... EACH ROW. WHEN ALL ROWS HAVE BEEN CONSIDERED THERE IS 

C... ONE VECTOR WHICH IS NORMAL (OKIHOGONAU TO ALL ROWS 

C... and CONTAINS THE SOLUTION VECTOR. NO MATRIX INVERSION 

C... IS INVOLVED and ONLY ONE ROW Of THE COEFFICIENT MATRIX 

C... IS REQUIRED AT A TIME. AND ONCE OPERATED ON CAN RE 

C... OVERWRITTEN. IN ADDITION. TWO AUXILIARY ROW VECTORS 

C... ARE NEEDED FOR TEMPORARY STORAGE OF INTERMEDIATE VECTOR 

C... VECTOR DATA, 

C . . . 

c . . « 


TO SOLVE THE LINEAR SYSTEM OF BOUNDARY CONDITION 
equations by PURCELL^S VECTOR ORTHOGONAL I ZAT ION 

method. 

CALLING SEQUENCER 

ITOTAL = TOTAL NUMBER OF HORSESHOE VORTICES. 

NXl = ITOTAL ♦ 1 

EW = ROW OF NORMALWASH INFLUENCE COEFFICIENT MATpIX. 

COMMONS 

ALOC. 

CALLING SEQUENCES 

VORl = auxiliary computational Row vector. 

VORK = extended SOLUTION ROW VECTOR. 

VORL = AUXILIARY CUMPUT A f TONAL ROW VECTOR. 

COMMONS 

GAMMA. 
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CONTROL^VRLX.WASH 9/14/ 

c • • • 

SOBROUriNt «|ASH (X, Y, Z, St T* H2t Ut V, Wt AAt A^t TEt rTt MN 

C . • . 

C... PURPOSE TO compute 1 HE THREE VELOriTY COMPONENTS INDUCED AT A 

Cf GIVEN POINT BY A GENCRALl/ED HORSESHOE VORTEX OF UNIT 

c... strength. 


CALLING SE0UENCE< 

At Y* Z = 0RTH0G0^JAL CAPfESlAK COORDINATES OF 

receiving (FIELD OR CONTROL) POINT MEASURED 
IN A REFERENCE FRAME Cf-NTEREO AT THE MfDPOlNT 
OF THE transverse VOwTpX SEGMENT (hORSrSHOE 
VORTEX CENTROID) + Tht X-AXIS IS PARALLrL TO 

The X-axis i)f The mastfr (C(jtsiF igurat lo u 

COORDINATE system* THE Y-AXIS IS NOt<iMA| TO 

the X-AXIs R<JT lies in the plane determined 

HY THE X-AXIS I T SELF AND THE TRaNSvERSf I EG 
OF THE HORSFSHOFt AKu THE Z-AXIS IS NO^/MAL 
TO SUCH PLANE. 

S = HORSESHOE VORTEX SFMISPAN. 

T = tangent of SWFEP AN(il.F OF TRANSVERSE I EG* POSITIVE 
FOR SWEEPBACK. 

8? = COMPRESS/ MIL I T Y FACTOR (= MACH - 1.0). 

AAt AM = direction ANGLES (ANCLES OF FLOATATION) OF 
FREE (WAKE) TRAlLiNCi LEGS. 

TE = TANGENT OF TRAILING ruGE SWEEP ANGLE. 

CT = AVERAGE LENGTH OF HOUND TRAILING LEGS OF HUpSFSHOE 
(DISTANCE BETWEEN VORTEX CEnTRUIO AND TRAILING 
EDGE MEASURED ALONG X-AXIS>, 

MM = FLOATING WAKE COMPUTATION F|AG. 

COMMONS 


C.. .OUTPUT 

C . • . 


CALLING SEOUENCE< 


= ORTHOGONAL VELOCITY COMPONENTS INi;>UCED BY 
generalized HORSESHOE \/ORTEX OF UNIT 
CIRCULATION INTFNSlTr. 


c. ..subroutines 

Ct* -CALLED NONE. 

C « • • 

C. • .0 1 SCUSS 1 ON The I GEFiERALIZE(J( horseshoe VORTfX ELEMENT CONSISTS OF 
€••• FIVE LFC»S OR SEGMENTS* OF WHICH THREE ARE I HOUND f AND 

1 WO ARE [FREE/ OR T FLOAT INGI. I HE HOUND LF(jS ARL THE 
C..* SKEWED* OR SWEPT* TRANSVERSE StGMENT AND THE TWO 

C*** TRAILINTi* or CHORDWISE* FJLAMEimTS EXTENDING FROM I hE 

C*t* ENDS OF The transverse LEG TO THF TRAILING EDGE* THE 

FLOATiNfi TRAILING LEGS ARF. THE SfMI-INFINITE LINfS 
C*** WHICH START AT T H£ TRAILING EDGE AND CONTINUE TO 

uownstrfam Infinity according io a prescribed oi-^lcmon 
oetermimed ry the floatahon angi.es aa and am. These 
L-** floating IRAILING legs const ITuTf The continuation 

IN The wake of the round trailing segments. 

c-*- AT SUPERSONIC MACH NUMBERS THE VfLOCITY INDUCED hY A 

C*** DISCRETE HORSESHOE VORTEX BECOMES VERY LARGE IN THE 

C*** VICINITY OF THE MACH AFfCONES oFnERATeD HY THE Sktwf^O 
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74/74 OPT«l 


TTN 4*b^4IOA 


10/04/76 lO.lS.l ^ 


LEO OF THE HORSESHOE. AT THE CHARACTERISTIC ENVELOPE VORLAX ‘564H 
ITSELFt The induced VELOCITY VANISHES DUE TO THE FINITE VORLAX 3649 
part concert, this singular behavior OCCURS ONLY FOR VORLAX 3670 
FIELD POINTS OFF THE PLANE OF Thf HORSESHOE. TO AVOID VORLAX 3671 

THIS numerical singularity the characteristic SURFACES VORLAX 367? 
(MACH CONES) ARE TAKEN TO BE 6IVFN BY VORLAX 3673 
(X - xn ••a « B2 ♦<(Y - YI) ♦ U -71) ♦♦?) /CUTOFF VORLAX 3674 
NMERE CUTOFF IS A NUMERICAL CONSTANT WHOSE VALUE IS VORLAX 367S 
LESS THAN* HUT CLOSE TO* 1.0 . VORLAX 3676 

VORLAX 3677 
VORLAX 3678 
VORLAX 367R 
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74/7A OPT=»l FTN A.5 + 410 a 

10/04/76 

10. IS 

control^vrlx 

•ZNORM 9/14/76 

VORLAX 

3B94 

c • « • 


VORLAX 

3895 

SOBROUTINE ZNORH (1» 1S« RTOP) 

VORLAX 

3896 

c • • • 


VORLAX 

3897 

C...Pt)RPOSE 

TO integrate the chorowise slope distribution in order 

VORLAX 

3898 

C • • • 

TO obtain the surface camber (OR RARP) ORDINATES. 

VORLAX 

3899 

c 


VORLAX 

3900 

C... INPUT 

CALLING sequences 

VORLAX 

3901 

C • • • 

I * ELEMENT INDEX (TOTAL COUNT). 

VORLAX 

390^ 

c • • • 

IS » STRIP INDEX (SPANWISF COLNTj. 

VORLAX 

390 3 

c • • • 

RNMAX s NUMBER OF HORSESHOE VORTICES In A GIVEN STRIP 

VORLAX 

3904 

c • • • 

(CHOROWISe ROW). 

VORLAX 

390S 

c • • • 

COMMONS 

vorlaa 

3906 

c • • • 

LAX* SLF» SLOPE. 

VORLAX 

1907 

c 


VORLAX 

3908 

C.* .OUTPUT 

CALLING sequences 

VORLAX 

39U9 

c * . . 

NONE. 

VORLAX 

191 0 

c * . . 

COMMON < 

VORLAX 

3911 

c • . • 

ZC. 

VORLAX 

3912 

c * . * 


VORLAX 

)91 3 

C... SUBROUTINES 

VORLAX 

3914 

C...CALLEO 

NONE. 

VORLAX 

3915 

C * • * 


VORLAX 

3916 

C. . .DISCUSSION THIS SUBROUTINE COMPUTES THE CaMrER NORMAL OWDINaIES 

VORLAX 

3917 

C . . . 

BY TRAPEZOIDAL INTEGRATION OF 1 HP SURFACE SLOPE OlS 

VORLAX 

391 8 

c . . . 

TRI0UTION ALONG THE CHORD, THIS SUBROUTINE IS ON| Y 

VORLAX 

3919 

c . . . 

called (HY SUBROUTINE PRINT) WHEN THE DESIGN PROrESS 

VORLAX 

3920 

c . * . 

IS INVOKED (INVERS = 1). IN SLCH A CASE^ THE LOAo niS 

VORLAX 

3921 

c . . . 

TRIBUTION HAS BEEN INt»UT* AND IHF SURFACE WARP NfEDED 

VORLAX 

3922 

c . . . 

TO ACHIEVE IT IS THE DESIRED OUTPUT. THE CAMBERLTNF IS 

VORL AX 

3923 

c . . . 

DEFINED AS BEING ZERO AT THE LtAOlNG EDGE AND IS 

VORLAX 

392 4 

c . . . 

expressed AS DECIMAL FRACTION OF THE LOCAL CHORD 

VORLAX 

3925 

c . . . 

length, the computed camofrlinf represents The Total 

VORLAX 

392^ 

c . . . 

SURFACE WARP* I.E.* IT INCLUDES CAMBER AND TWIST. 

VORLAX 

392 T 

c . . . 

SUBROUTINE ZNORM IS CALLED ONCE PER EACH CHOROWlsE 

VORLAX 

392 6 

c . . . 

STRIP OR ROW OF HORSESHOE VORT ICPS. 

VORLAX 

3929 

c 


VORLAX 

39 30 

c 


VORLAX 

39 3 1 

c 


VORLAX 

39 32 
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STATISTICS 

PROGRAM LENGTH 2750B 1512 

CM labeled common length 7711 IB 32329 
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ALPHA = 12' 


OUTPUT DATA 


DELETED HERE 



APPENDIX C 


WAVE DRAG TO VORLAX INPUT 
CONVERSION PROGRAM 



SUMMRY 


The purpose of a program called WDTVOR, developed to convert the Wave Drag 
input geometry into the VORLAX input geometry description, is to save time, 
improve accuracy, and reduce human drudgery when configurations for which the 
geometry was first digitized in the Wave Drag format are also to be analyzed 
on the VORLAX program. 

The present version of WDTVOR contains the option to convert fuselages to 
flat plates having the correct planform area or to a simulation having hexa- 
gonal cross sections. All wings and planer surfaces are converted to zero 
thickness panels although the wing camber effects are preserved. Engine pods 
are converted as curved panels approximated by hexagons. 


INTRODUCTION 


In the analysis of an airplane design, what is basically the same geometric 
body is described by several different geometric models, each of which is unique 
to the discipline for which it is designed; yet at the same time there exist 
certain elements common to all models. For example, the NASA Wave Drag format 
(Reference 12) emphasizes the enclosed volume of the aircraft and correct spatial 
relationship of components while the VORLAX method of aerodynamic analysis 
(unified subsonic and supersonic) uses a paneling procedure similar to but 
different in detail from the NASA Wave Drag format in order to reduce computa- 
tion time. 

The program described herein, called WDTVOR, is for direct transformation 
from WAVE DRAG to VORLAX input format, A comparison program, called VORTWD, 
transforms data from VORLAX to WAVE DRAG formats and is documented in 
Appendix D. 

Presently the simulation form into which Wave Drag data is converted is 
limited to that most commonly used in VORLAX analysis as follows: 

GENERAL 


Cosine law spacing of streamwise vortices and equal interval spacing of 
lateral vortices is assumed. 
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A single Mach number from the first wave drag case card is picked up and 
output in the VORLAX dataset. Three angles-of-attack, 0, 5, and 10 are also 
included. 

Symmetric P'^light is assumed. 

The value of CBAR is arbitrarily set to 1,0. The moment center coordinates, 
XBAR and ZBAR, are set to zero. 

All major planar panels are set up for three spanwise vortices (NV0R = 3) 
and 10 chordwise vortices (RNCV = 10). Full leading edge suction is assumed 
on each panel (SPC =1.0) except for the fuselage panel. 

WINGS 


• Wing camber is preserved in the transformation. 

• Wing thickness is not preserved. Converted wing panels are zero 
thickness . 


FUSELAGE 

• Both circular and noncircular cross section fuselages may be converted 
to a single zero thickness panel of trapezoidal shape. The planform 
area of the zero thickness fuselaTge is made equivalent to the planform 
area of the Wave Drag simulation and the vertical location of the 
fuselage is set equal to that of the leading edge of the wing root 
chord so as to effectively "seal" the wing-fuselage juncture in the 
VORLAX simulation. 

• Optionally fuselages which are circular in the Wave Drag simulation 
may be transformed to a hexagonal cross-section simulation in VORLAX. 
With this option transverse vortices are located on the actual body 
surface. Also a wing- fuselage seal is assured by holding the fuselage 
diameter constant in the region from leading to trailing edge of the 
first wing chord and by assuming that one vortex of the reference 
P®lySon lies in the plane of the first wing chord. 

® 1^ the Wave Drag simulation was a cembered circular one, the fuselage 

camber is preserved in the VORLAX simulation, 

• If the fuselage was of abritrary cross section in Wave Drag no camber 
information is transferred to VORLAX. At present it is possible to 
convert arbitrary fuselages into only the zero thickness planform 
simulation in VORLAX. 


NACELLES 
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• Circular nacelles are transformed as closed hexagonal panels using the 
option that transverse vortices are located on the actual body surface. 



FINS 


• Each fin is treated as a single panel of zero thickness. 
CANARDS 


• Each canard is treated as a single panel of zero thickness. 

After a Wave Drag data set has been converted to VORLAX format it may be 
validated by converting it back into Wave Drag format and a configuration plot 
made. This plot will, of course, be different from a plot of the original 
Wave Drag dataset by virtue of the zero thickness panel representations. 


NOMENCLATURE 


CPU Central Processor Unit. 

FUSTYP Flag for fuselage type. Use 0 for single trapezoidal panel 
fuselage simulation. Use 1 for hexagonal (curved) panel 
fuselage simulation. 

VORLAX A Unified Non-Planar Vortex Lattice Method for Subsonic and 
Supersonic Flow Described in this report. 

WDTVOR Wave Drag to VORLAX data conversion Program described in this 
appendix. 

VORTWD VORLAX to Wave Drag data conversion program described in 
Appendix D. 


PROCEDURE 


1. Prepare a dataset in WAVE DRAG format if it does not already exist. See 
Reference 12 for the WAVE DRAG input description. Be sure at least one 
case card is included with the dataset. It is necessary to instruct the 
conversion program as to whether the VORLAX simulation of the fuselage is 
to be a single trapezoidal panel or a curved (hexagonal) panel. This is 
accomplished by appending to the front of the WAVE DRAG dataset a card in 
namelist format specifying either FUSTYP = 0 for a single trapezoidal or 
FUSTYP = 1 for curved (hexagonal) simulation. Figure C-1 is a listing of 
an example WAVE DRAG dataset showing how the namelist input is placed at 
its beginning. Figure C-2 is a graphic representation of the Figure C-1 
dataset obtained from the Configuration Plot Program of Reference 12. 
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2 . 


Submit the dataset with the WDTVOR program. A compilation listing of the 
program is included at the end of this appendix. 


3. Examine the output from the above submittal. Output consists of a listing 
of the input dataset in card image form (see Figure C-1) followed by a list- 
ing of the output dataset in the VORLAX format. Figures C-3 and C-8 are 
sample datasets as output from WDTVOR program. Figure C-3 was based on the 
ciirved (hexagonal) panel VORLAX fuselage simulation while Figure C-8 was 
derived from the single trapdzoidal panel fuselage. 

1+. Because of the simplification that takes place in WDTVOR, e.g., zero 

thickness panels, it may be desirable to generate a graphic representation 
of the VORLAX dataset. This is accomplished by using the VORLAX to WAVE 
DRAG conversion program described in Appendix B followed by the Configura- 
tion Plot Program described in Reference 12. Figure C-4 is a sample of a 
simplified dataset after transformation back into WAVE DRAG format. Figure 
C-5 is a plot of the Figure C-4 dataset. Figure C-6 is a listing of the 
WAVE DRAG format for the single trapezoidal panel fuselage simulation and 
Figure C-7 is the corresponding plot. 

5. Generally, a certain amount of editing is required to prepare a VORLAX 
dataset generated by transformation via the WDTVOR program for an actual 
run on VORLAX. Figure C-8 is a listing of a converted dataset as output 
from the conversion program and Figure C-9 is the same dataset after it 
has been edited preparatory to a VORLAX run. 

An explanation of the notes of Figure C-9 follows: 

NOTE 1 - The value of the over relaxation parameter, REXPAR was 
changed from -0.02 to 0.0 implying that the program is 
to compute internally the optimum over relaxation value. 

NOTE 2 - The number of angles of attack, NALPHA, was reduced from 
3 to 1 and the single ALPHA selected was 5.0 degrees. 

NOTE 3 - The number of major panels, NPAN, was changed from 12 to 13. 

The pitching moment reference length, CBAR, was changed 
from 1.0 to 1022.28. The XBAR and ZBAR of the moment 
reference point were changed from zero to 1910. and - 800 .i+, 
respectively. 

NOTE k - The number of spanwise vortices, NVOR, and the number of 

chordwise vortices, RNCV, were adjusted to provide a better 
subpaneling representation. 

NOTE 5 - Wing Panel Number 5 was completely eliminated. It was actually 
a null panel having little or no span which was included in the 
original WAVE DRAG simulation for another purpose but which 
would have caused errors in VORLAX if it has been retained. 
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NOTE 6 - A vertical panel vas added to simulate a nacelle pylon which 
was not included in the original WAVE DRAG simulation. 

NOTE 7 - The spanwise location of the nacelle was moved outboard from 
2lU. to 233.9^ to coincide with a major panel break in the 
wing. Alternatives to this procedure would be to create an 
extra major panel having an edge at the 21^4-. span location 
or making sure that one of the wing vortices resulting from 
subpaneling coincides precisely with the engine pylon span- 
wise location. The above precautions must be taken whenever 
major panels intersect. 

NOTE 8 - The variable area nacelles were changed to fixed area ones by 
setting the ZCI values to 100. This makes each nacelle cross- 
section eciual to that of the reference polygon, in this case 
the inlet area. This was a compromise made to simplify the 
nacelle pylon simulation. With constant area nacelles the 
pylon could be simulated by a single rectangular panel. 
Retention of the variable area nacelles would have required 
several panels to seal the pylon area between the wing and 
engine, 

NOTE 9 - The location of the nacelle was moved outboard from 20U. to 
233.9^5 otherwise, Note 7 applies. 

NOTE 10 - Fin Panel #2 which was on the aircraft centerline was 

completely eliminated. Since the desired VORLAX run was to 
be symmetric about the pitch plane the centerline fin was 
not required. Furthermore, its inclusion would have un- 
necessarily increased computing cost because its presence 
alone triggers off certain asymmetric calculations. 

NOTE 11 - The horizontal tail (labeled canard) root chord dimensions 
were changed to coincide with the side of the fuselage. 
Actually, only that portion of the horizontal tail outboard 
of the VORLAX fuselage simulation is included, 

NOTE 12 - A vertical panel was added to seal the gap between the 
plane of the horizontal tail (Z = -Uo.) and that of the 
fuselage (Z = 5^.96). 

6. It is a good idea to submit the edited dataset to the VORTWD program 
(Appendix D) to recheck the editing just done. This step is optional. 

It probably wouldn’t be done if only minor editing was involved in Step 5. 

If major editing was done, such as the addition of pylons and sealer panels 
as in Figure 0-9, then it is recommended. The results of converting and 
plotting the dataset of Figure C-9 is shown in Figure C-10, Note the 
absence of centerline fin, the addition of engine pylons and horizontal 
tail sealer panels and the constant diameter nacelle simulation. 
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7. The edited dataset may then be submitted to the VORLAX program. A portion 
of the VORLAX output resulting from this submittal is presented in Figure 
C-11 . 
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Figure C-1.- Typical WAVE DRAG format dataset 
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Figure C-1.- Concluded 
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Figure C-2.- Graphic representation of WAVE DRAG dataset 



o 


o 

o 

o 

o 


vO 

OQ 

<XJ 


u 

^0 


o eg 

^ I 


o 

o 

o 

o 

• 

o 


a 

o 

o 


tTi 

OC 

r- 

o 

• 

0 

1 


o 

o 

o 

o 

• 

o 

un 


Oj 

o 

• 

0 

1 


g: 

(\i 

eg 

O 

• 

0 

1 


o 

o 

o 

o 

# 

o 


o 

04 

O 

• 

0 

1 


CO 

eg 


o 


o 


o 

o 


*~l C) 

• • •» 

ooo 

2T2CLO 

00 

>0 J-J » 
^ujai> 

o<>2eo 
•0<<2 
O •O.O.QC 

o •* 
ooar 
04220 


o 

o 

c 

o 

• 

o 




o 


lO 

04 

CO 


04 O 
I i 


O 0) O 

o fo 
O fO 
O a og 

• • • 

° T ? 


o 

•o 
o • 

o 


OO00 4*sf O 

OOeg^^f-iN^ 
• ••••• 

OOf^ 4^000 
040 ( I I 


egogO 

a 

OOU 

2eao 

•> 

U4Ui> 

22TO 

5<2 

ooa: 

220 

3C3t2 


O 4T GO 

O 04 

O 00 14'. 

o >c ^ 

• • • 

^ ° 


• ro mo 
c a 

OOU 
220.0 
uo 

o c -j j » 

O fO 0- uluj > 
o 4“ eg 220 

04 m <<2 


• • 

o o 
I 


o 

4 


aatt; 


oo^ftf^oor^ 
OO4‘O«"0f^ 
OO«*4Q0<M ^ 

• ••••• 

OOC<OiOi-i 
I I 


220 
H4^H>0 
:S3t2 


00 00 
cn e- 

ro 


o 

o 


44C 

a 

OOo 

220LO 


OO OOf^iT^CNJ 
OOroh-^^ 
OO Oco 
004 >^fOrg 

• ••••• 

ooo#^o<^ 

NO I I 


220 
:r:52 


•13 

O- 


o 

o 

o 
4 


lO 

04 

a 


o o 


0 

r- 

0* 

io 

uu 

0 

O'* 

en 

0 

0* 

>0 

Ui UJ> 

0 


cc 

0 

eg 

\T\ 

220 

0 

4n 

eg 

0 

n 

eg 

<<2 

0 

eg 

pH 

• 

• 

• 

ixo^ae 

• 

• 

• 

0 

0 

0 

»■ 

0 

0 

0 

fO 



oucr 

n 




OO 

ooO(nif\rn 

OOroogom 

OCOi-400 
040 I 


a 

ujO 

2rg 

•-40 

o • 

5? 


eg 


u 

> 


2 

oe 


O O'O 

o c^o 

OO no* 

• •OO'O^O 

CO • • • • 

O004T»O 

-404 

4C0 

P-4 CO 


ocroocoee4 0oo 
00^00^0" o<o 
OO®OO'C0*O 

o0 4cono^oo 

• •••••••• 

ooonon^joo 
I I I 04.n 
7' 43 


rn 


OOooegog#^GC<J' 
OOO*^ 7 

OOO'04NO'0 
oofoonno oco 

000r0o*--<00 0 

t I eg 


O O eg <.0 o 

OO ^ n <N o ^ o 
004eg>0'C'004 
OOrOrOegO 00*0 

•••♦••••# 

ooo^aq>p^^o 


o oooo 
o oooo 
OOO'OO^O 
• *0^40 

iT^O • • • • 

••■•4 4 »"*4 

lAiTk 

I I 


fO 


T (4r-l 
.nro 


OOn Oirrg 
OOegOg 
oo^^r^eg 
CO 040^ ^rn 

• ••••• 

OOOOQO 

^7 i T 


000*0)004000 
pOsO'OOfOOOO 
O0r0s0r*rgr-00 
OO*-^ ^04^4 00 

*ei o on 004 4 0 ^ e-4 

OC I I I >c 
I I 


n 


OOO04 00 >0 o o o 

poon 04000 

Ooo^rgew ooO O 

O o eg sO «-H o »-<0 

• •••••••• 

«ooorgooar^»H 

OO I I Ofw 
I i 


00404rno4 
O O r- o 04 pH 
O O u rH 

oooo* r-in 

‘l^OOOOO 
00 I I 


2 

< 


4 

eg 

I 

wO 


O 

»T 

iTi 

* 

04 


oo<x)>c 4*gooo 

OC040-0000 
OOf^X0*Og34O 
OO'-HO' OO* o 

iTiOOpOQO*a r-* 
ou i I O^oo 

o ‘ ' O 

oo«-4-4-o ^oc ooS^^2®o?ic 

0003 ••••».... . . , , , , .O-^OO 4■OOOCJOO^0OvrOa3OOOOO^t0t^ 


g3 


eg 


fgpHr-4 

ogro 


i 

nn 


I 


00 hO 00*-hooO00 


i g3co 
n4 


I I 


S^oo o ^ 5oao o 2 o f- 

<>00 O O i?5o^O O 2 I^O'O'O o eg m 

oinooooo»HO'04A0*OQoSr'i?J,'^i^£.? ^ ego^ao o eg ^ 

3r=°rv°T|r°°s°r-?2S"°“s-f>ss— 

-iOpHcnoeg #-i O^eg 

U iH pHrg 


eg 

n 

I 

0* 

o 

g3 


260 


I'igure C-3. Listing of output dataset in VORLAX format 
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Figure C- 3 .- Continued 
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Figure C-3.- Concluded 
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Figure C-^.^ Result of converting circular fuselage VORLAX dataset tack to WAVE DRAG format 
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Figure C-4.- Concluded 
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dataset back to WAVE DRAG format 



fSjCSlfVJ 


O OC OO OO OO CO o 

or tt:or: ortf oca 

OCOOZOOZOCZOOZoBZO 
Q Q O ? 2T r Z r 2 2: < Z 

o a. O U.I*4 M LL u. IM ^ LL l-< < O < 

a X O. LL XU.U. K U.U. XU. IL >r Ul OKU 


oo 

« « 

CO 


oo 

• • 

oo 


0 

0 


CO 



00 • 

• 

• 

• 

• 

K-K 

• tfO 


•H 

in 


QCor 

co<^ 

0 


r«*i 


00 


m 

i-i 

CM 


cm 






0 (M 






• • 


0 

0 

0 

0 

0 (M 

00 • 


• 

• 

• 


• 

CM 

00 

in 

0 


ooir. 

r- 

fvi 

in 



1 

I 

1 

1 

1 



0 

0 

f- 



00 • 

• 


• 

• 


• #00 

m 

in 

in 

>r 


OOf\» 

C>l 

Nt 

>0 

r- 


<Vi 




f— < 


0 

0 


>r 

0 


00 • 

• 

• 

• 

■ 


« 

rH 

xCi 


-r 


000 

0 

0 

r- 

m 


0 

0 


*-» 

m 


ri 

(M 

rj 

fO 

fO 


0 

0 

in 

OD 

00 


00 • 

• 

• 

• 

• 


• *fO 



m 

in 


oo<> 

0 

0 

m 

m 


fO 

m 

t— 1 

<M 

CM 







m 

0000 


0 

0 

0 

(M 

000 • 

• 

• 

• 

• 

• 

• • 

- 4 - 

0 

0 

0 

o 


tn 

0 

>r 


cn 

oacn 1 

1 

f-H 

1 

1 


r-«fO 


1 



• 

1 





4 - 

^'to<>ooo^ooooor^cor^oo (nj 

0^0 • • 

• • • 

• • • 

• • • 

• * • 

• 1 

• • •<noof 0 OoiriOO 4 AOO»rkOO <m 

rosOCvJogO 

(MO 

^0 

M )0 

'OO 

M) 

mCMr-l 

(M^ 

<0^ 



• 

CM 





M) 






CM 

0 0000000000>t00't 00 

000 • • 

• • • 

• • • 

• • • 

• • • 

• 

• • • st 00 r~i 00000 *-* 0 o*-^oo 

•^0 (MO 

0 

CD 




0 fOO 

0 




>IN| 

0 (M 

(M 

CM 

fO 

<n 


(\i 





XX 


267 


Figure C-6.- Concluded 
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Figure C-7.- Plot of figure C-6 dataset 


o 


o 

o 


o 

o 


o 

t 

o 



0 

<M 


0 


W . 

0 

u\ 



0 

h 3 

0 


0 

>0 

0 

0 

0 

00 

0 

00 

(M 


0 

CM 

pH 


0 

CO 


0 

>o 

r- 

0 

p- 

CM 


0 


r») 


0 

CO 

pH 

0 

PM 

0 

0 

0 

0 


0 

0 

rvj 


• 

• 

• 

• 

• 

ft 

ft 

ft 

ft 


ft 

ft 

ft 


0 

tM 

rH 

0 

P-x 

0 

0 

a 

0 


u 

0 

0 



1 

1 

0 

1 

1 

iT 

1 

1 


tn 

1 



0 

0 

vn 

0 

yT' 

CO 

0 

CO 

o: 


0 

O' 

mA 


0 


(\! 

0 


0 

0 

r i 



0 


CM 


0 


00 

0 

CO 


0 

UA 

iO 


0 

n 

>> 


0 

>r 

sO -J 

0 

>0 

fH ^ 

0 

•-H 

F-* 


0 


pH 

p-ipH 0 

« 

• 

• M rvi 

• 

• 


ft 

ft 

ft 


ft 

ft 

ft' 

Q_ 

0 

rsi 

0 CL 

0 

0 

0 Q- 

0 

0 



'O 

0 

0 

« • •» 

nT 

1 

1 • • » 

-t 

1 

ft ft *• 

4 



ft ft •> 

>r 



CJQU 



OCJO 



cuo 




aaco 






2 ^ 0.0 


220.0 



220.0 



00 






CO 




UO 





0 

n 

0 -J-J •* 

0 

0 

pM 00 » 

0 

1 ^ 


-JO ► 

0 


P. 

0 ^ uJUJ> 

0 


n ^uj> 

0 


r* wuuj> 

0 


0 

LU»^> 

0 

«0 

CO 

oa-zeo 

0 

n 


0 


<M 22 U 

0 



220 

0 

\r 

<M 

•o<<z 

o 

O' 

fM <<2 

tJ 

fvj 

^<2 

0 

0 . 


<<2 

0 

:M 

pH 

0 *a.a.i^ 

t 

• 

• Q. CLaC 

• 

• 

• OlXX 

ft 

ft 

ft 

IXOX 

ft 

ft 

ft 

0 •• 

0 

pH 

0 ► 

0 

0 

0 ^ 

0 

0 

C 

•» 

0 

0 

0 

T^COx 

m 

1 

1 oox 


1 

OOX 

(n 



oox 

f') 






2^0 



22a 




220 














jc:* 2 



31X2: 



Jt32 




3 JC2 





00 

C' ■*T 

oo>^h»oor^ 

00 » rj 


00 

^ACMfrjOn 


0 0 CO xt Q 

oo^oc^r- 





0 

0 0 U vC go 

0 0*^0; 

i-<0 

OOMCOt'X 


0 0 0 

•0 

OO^pHi-i>r 


oo*4‘'Cr^ivi 


ooncviocn 

0 » 

t • 

• 

• • • 

• 

• ft • 

ft ft 

ft 

ft ft 

ft ft ft 


ft ft 

ft 

ft ft ft 


om 

rsic I I I 


aoono«-i 
I Vi o I i 


^^0 I i 


ooo-<oo 

‘MO I 


O 

OjO 

t~»o 

o • 

2:0 
^ I 


Cs 

o 

> 


a: 

or 

OrH 


o 

o ;^o 
00 

• •OC'^O 

00 


oonccaoMOOj 

OOO'Or^ooch 

ooooac)>o^^o 

oO'tomOvj'OO 


o o (Si tNi »-« tn 

OOO'U 4‘ ViO< 
OOmoroff*.oano 


ooM«-i,n 4j a o 
COa rn>-s40<? O 
co-r-'j-a^^cr j 
OOi'OrOogOvOU^'O 


00(M0u^f*- 

CO 


^ OcOJiTiO s-^OO^^'iOo^iTiOv^ OOOC^O»-<O^Q C O ^-*5 >-<0 O yj' piO OOOOOO 

'tM »^C' I I I (MU^ p-iU' I I iTuM ( i < I 

O >0 u*> ir tr, 

n n. r'i (<. 

^ 1 -^ ^ « 

o 0000 oor^ococsdooo ooocmoowdooo oooo<^fViOoo co>trNjrnr>i 

o 0000 ooso^oroooo ooof^o 4 ooo ooostr^ocoo oor^O(\iiH 

ooo^h*o oom^or^csih - 00 oo^<Mr^ooo^o ooi^o f^oo^to ooh-o<-4rH 
• •oa»to ooi^»crsj>o>t 00 oor\i^p-«o*'0^o oo»-<0"0f^‘-^r^o ooor^fHfo 


a^O 


ITOOfOOrM'ta r-J U^OOfMOOCM*^r-t inOOOOO(^>F^ LTOOOOO 
OD I I I » I I T Ti^o a T i 

t i I ( I I II 

0000 


00 T<>coina 0:0 O0ooir,p-«^r>oj0 
oocMs^mcoa O'O O0f0oovni-ic>a o ooir«-iir.r-i 
oOif r*^w'Ot^o r »r^o oOi-iu‘,f\.cr 


2 O 0000 OOv^OC^O'Oao 

< iC\ 0000 OQ<Mf<'i<MNOOO^O 

-J tnOO •vf^O OOi-H>t»rr^^OO _ 

a. • • •Oi^a^OOOOO^Oi«> ^OOOOOf^O^^OOOOOO^O i'Oa.OOOOO >f 0 <N 
c^^coo: 

>riOfsimr-co^oocgooi^ccQO.pHOooopcQsOOO-HOoooo 
I I 1 I I T Ocori T I 

<Mn rO>t 


(M 

I 


'O 00 r^O p < 

f*^ >C iM »H I 

*C (VI 


JsOO 

^00 

^00 


(M 
n*i 

I 

o 

sO 

*-* O'st 

-JOiHfOOrg #H 

u 


o 

o 

o 


4 

cn 

O' 


u^Oo^O 

4 -OChO 


o 

o 

c 


iiT 

4000^0 

0 


O' coo 

0 

p- 

vP 

ino^O'Q 

0 


rsiO'OQ 

0 

(M 

U'^ 

u>^OvO 

0 

fcA 


0 

(M 


OknoOOOO«-iO'^^h“OOOOO^Or‘’>h-mOOOOO<<^Orsir l^"000 00<MO^N 
'4*<MnoOOOrnCF-*r400rOOOOOi-iOOOOrf^mOOOOOOO*r»U^t^OOOOOOO 

^4^x0 >0 I 1^^ nO I 1^0 o T Tofvi o T 

'4‘^• r^C' O^CM 

F-IfNj 


269 


Figure C-8.- Converted dataset before editing for VORLAX run 
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Figure C-8,- Continued 
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Figure C-8.- Concluded 



CL 1607-32 (-24 PLANFORM) VC 121 EN(? 1 NE NOTE 1 


o 


Csj 

O' 

0 

O' 

in 

0 

m 


>C 

0 

0 

0 

00 

0 

QO 

<M 

CO 

0 

0 

'O 


0 

n- 

fsi 

09 


0 


0 

0 

0 

0 

• 

• 

• 

• 

• 

• 

• 

• 

fvj 


0 


9 

0 

9 

9 

1 

1 

m 

1 

T 

in 

T 

T 


^ c 
o. 

o o u 

2 2 a O 
</> 

%0 -J -J 

a Ui Ui > 

o a r r o 

UJ • O < < 2 

u. O • Cl O. Ct 

O O 

7 O O a: 

<M 2 Z O 

V ^ M IM > o 

V- * » 2 

• 

o o 

• o 

O 09 

I 


0 

0 

m 


0 

in 

00 


0 

00 

CO 

0 

K. 



0 

rvi 

fO 


0 

cn 


8 

K 

09 


S 

00 

in 


0 

in 


0 

4 

nO 

•J 

0 



•J 

0 



• 

• 

• 

M rg 0 

• 

• 

• 

fo 0 

• 

• 

• 

0 

cu 

1 

0 

CL 

0 

0 

0 

(L 

0 

0 

0 

4 

• 


4 

i 


t • ^ 

4 





OOO 




ooo 





2 Z CL 


cn 

0 

ul 

0 0 

0 

0 


0 

(/) 
•J ^ 

0 



h- 

fO 

UJ Ui > 

0 

fO 

1 ^ 

UI 

Ui > 

0 


'O 

m 

4 

2 2 0 

0 

4 

m 

2 

2 0 

0 

<M 

tn 

0 

ni 

< < Z 

0 

<M 

fO 

< 

< 2 

0 

fO 


• 

• 

a CL O' 

• 

• 

• 

a 

(L OC 

• 

• 

• 

1 

0- 

» 

0 0 oc 
2 2 0 
H* >-< > 0 
:s or 2 

0 

fO 

? 

0 

WING 

0 OC 

2 0 

M > 
!]B 2 

0 

m 

0 

0 

0 


O O O' 
O O 00 
O O in 
O O M 

• • • 

O O r4 

rvj o I 


O iO 00 
o 4 - o m 

O ^ 00 r-( 
O r-l ^ 4* 

• • • • 

o o ro o 
o I 


O 00 in <M 

o <n O' ^ 
o ^ o fO 
o ^ 'O m <Vi 

• • • • • 

o o •-< o ^ 
O I I 



^ O' o 

• O' o 

o o 

t ^ a O' o 

O O' * * * 

•H oc in o 

— ' 4 * (Ni 
O' 


oofncooon^ow 

8 0 0 o r- O o O' 
o 00 o o ^ o 

OO^O'mOO'OO 

ooorncfninoo 
o I I I M in 
o o 


OOoofMCsJ^OCO' 

oor^O'O'foO'O' 

OOO^'tCViO'O 

oo<nomcooooo 

OOO(<>Of- 4 O 0 'O 
^ o I t «n esi 
in 


o o <v on ^ 

o o O' fo nj o 

o o ^ ^ < 

O O fO fO fvi o 

• ••••• 

o o o o o 

•-to I T 


09000 oor^flooniooo ooofvjoooooo oooo^^^ 
<MOOo*^oO'^'OOfnooO'^ooomo^tooO'^ooo>tf'-o 
o •o^•OLiJOO<n^^•N^“OOLuOor-^^J^•ooo^oOliJOO^•op^p 

• nj 0 '^ 0 u- 00 f-i^ivi' 04 ’ 00 )_ 00 ni* 0 «-t^ 0 «- 40 )- 00 r-i 00 r^ 


IfVl^tCQ***** 

| 0 ' 4 ''t*-^ 7 in 00(<>0 

iip in oil 


• • • 

n^ O' 

' Tf 


• • o • 

h ^ »n 


yinoorNiOOO'h'»~i7inoooo 

^ I i m i 


or*-t oooo OOOO'O'O'OO'O ooo'o'oomo'floo oooomi-<^ 

in oooo oo(Mf 0 fsl^o<^o oonj^focoO'O'O oomooin*-! 

moo •' 4 '^c 00i-<^m^'4‘00 oom’4'«-^moroo oo*-iin«-im 

• • •or'- 0 'Ooooo>fo»-HO'tooooo»-<o^'toooooo 4 )O 4 ' 

fvjmooQ^** •*•*•••••♦••••••••• •••••••• 

LI_ >omm'Cio*"ioOfoo<Mfnh*>oOp-<oof^oof^<x)' 00 «^oooop 
>C ^•ll I I I 


I <\J i-» 

<N4 fO 


ooo o ^ moooo o m 4 'odo'O o ^ f'* 

O'OO o m ^toO'O o 4“ mO'O'O o m ni 

O'OO o O' inooo o m mo'O'O o m m 

omOoooo*>^o^mr-ooooo%oOfnr'-rnooooofnonj 

4 fsi 40 oo oc'io*-<<Mcomooco^oooDmrsiooc50poo 

4 >c I I 'C >cLJ I T 'O ouj 'Oil 


272 


Figure C-9.- Converted dataset after editing for VORLAX run 
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Figure C-9.- Continued 
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Figure C-9.- Continued 
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Figure C-9«- Concluded 
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Figure C-10.- Plot of dataset edited for submittal to VORLAX 
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Figure C-11.- Portion of output from submittal to VORLAX 
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APPENDIX D 


VORLAX TO WAVE DRAG INPUT CONVERSION 
PROGRAM (VORTWD) 




SUMMARY 


A program developed to convert the VORLAX input geometry description into 
the Wave Drag input geometry description has two purposes: l) to permit plot- 

ting of the configuration geometry in wire frame form as a check on input 
errors; 2) to save time and reduce human drudgery when configurations for 
which the geometry was first digitized in the VORLAX format is also to he 
analyzed for wave drag. 

While the present version of VORTWD does not convert all VORLAX input 
options, it does handle the most common ones. It is recommended that all newly 
created VORLAX data sets be converted and plotted to validate the input geometry. 


INTRODUCTION 


In the analysis of an airplane design, what is basically the same 
geometric body is described by several different geometric models, each of 
which is unique to the discipline for which it is designed; yet at the same 
time there exist certain elements common to all models. For example, the NASA 
Wave Drag format (Reference 12) emphasizes the enclosed volume of the air- 
craft and correct spatial relationship of components while the VORLAX method 
of aerodynamic analysis uses a paneling procedure similar to but different 
in detail from the Wave Drag format. 


The program described herein, called VORTWD is for direct transformation 
from VORLAX to Wave Drag input format. To accomplish this, certain data mus 
be added to the VORLAX dataset to give the program sufficient smarts as to 
which VORLAX panels are associated with which V/ave Drag components. 


A companion program called WDTVOR converts data from the Wave Drag to 
VORLAX formats and is described in Appendix C. 

Presently the VORTWD program is limited as to the types of VORLAX data 
sets it can handle as follows: 


VfINGS 

• Wings may have camber but not thickness. 

• All wing panels must have the same number (NAP) and values (XAP) of 
percent chord stations for the camber definition. Flat wings without 
camber can also be handled. 
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• All wing panels must span the distance from leading to trailing edge 
of the wing. That is, separate panels representing leading or trail- 
ing edge flaps are not properly treated at present. 

• Twist applied to a panel via AINCl and AINC2 in VORLAX is ignored in 
the transformation. 

FUSELAGE 

• VORLAX fuselage simulations consisting of either a single trapezoidal 
panel of zero thickness or a curved panel having polyagonal cross 
sections may be converted. 

NACELLES 

• Circular nacelles represented by curved panels in VORLAX can be 
converted. Nine is the maximiam number of nacelles . Nacelles 
treated as vertical flat plates are better treated as fins . 

FINS 

• Zero thickness panels without camber are treated; one panel per fin 
maximiim. Six is the maximum number of fins and fin panels. 

CANARDS (Horizontal Tail) 

• Zero thickness panels without camber are treated; one panel per 
canard maximum. Two is the minimum number of canards and canard 
panels . 

Two plot control cards compatible with the Configuration Plot Program 

(Reference l) are automatically added to the converted dataset. 


NOMENCLATURE 


Symbol Description 

NWNGP Total number of VORLAX panels making up the wing description. If 
wing is cambered, then all panels must have the same number of 
percent chords (NAP) and the same values of percent chords (XAF) 
in the camber description. This is a Wave Drag requirement but 
not a VORLAX requirement. Input as an integer, no decimal point. 
Maximum value is 20. 
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Symbol 


Description 

NWING VORLAX panel numbers making up the wing. VORLAX panels are 

numbered according to their input sequence, not by any alpha-- 
numeric description appearing to the r_,-,ht of the VORLAX data 
cards. Wing panel numbers need not be consecutive but they must 
be listed in a sequence proceeding from root to tip. Input as 

integers separated by commas. The maximum number of values is 

20 and the number listed must be equal to the I'lWNGP value. 

NFUSP The total number of panels making up the fuselage. The fuselage 

is assumed to consist of only one major fuselage segment in terms 
of Wave Drag. Input one as an integer. 

NFUSE The VORLAX panel number describing the fuselage. 

NP Number of curved panels describing pods or engine nacelles . 

Maximum value is 9» 

NPPAN List of VORLAX panel numbers making up the pod or nacelle descrip- 

tions, One panel per pod. Presently, such panels should be 
curved panels, i.e., PDL>360. If a nacelle is being simulated 
as a vertical flat plate or plates, then it would be more appro- 
priately treated as a fin or fins for Wave Drag purposes . 

NF Number of fins to be described. Symmetric fins need be described 

only once. Maximum value is 6, 

NFPAN List of VORLAX panel numbers making up the fin descriptions. 

Maximum number is 6 and number listed must be equal to the NF 
value. That is, one panel per fin. 

NCAN Number of canards (horizontal tails) to be described. Maximum 

value is 2. 

NCANP List of VORLAX panel numbers making up the canard (horizontal tail) 

descriptions . Maximum number is 2 and number listed must be equal 
to NCAN. That is one panel per canard. 

PLOTSZ The maximum dimension in inches of the plot to be made from the 
plot cards appended to the converted dataset. Same definition 
as PLOTSZ in the Configuration Plot program (Reference l) . 


PROCEDURE 


1. Create a VORLAX dataset if it does not exist already. 
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2, 


Add smarts cards to the front of the VORLAX dataset. A typical dataset 
with smarts cards included is shown in Figure D-1. Note that all such 
cards are input in namelist format in which nothing may be entered in 
column one; only columns 2 through 80 may be used. The first namelist 
card must contain &INPUT and the last card &END. Input quantities are 
input as VARIABLE NAME = VALUE, as shown in Figure D— 1. Variables may 
be named in any order. Definitions for the various smarts variables 
are found in the nomenclature section. 

3- Submit smartened dataset to the VORTWD conversion program. Figure D-2 
shows the output resulting from conversion to Wave Drag format. 

Validate the converted dataset by submitting it to the Configuration Plot 
program (Reference l) . Figure D-3 and D-1* shows the plotted results of 
the Figure D-2 dataset. 

5. If the configuration is to be analyzed for Wave Drag; the converted data- 
set must be edited to supply thickness distributions for the wing and 
tail surfaces and if the VORLAX fuselage simulation was the single 
trapezoidal kind, then it must be given volume. 
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INPUT DATA LISTING FOLLOWS 
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Figure D-1.- VORLAX dataset with smarts cards included 
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Figure D-1.- Continued 
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Figure D-1.- Concluded 
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Figure D-2.“ WAVE DRAG format dataset resulting from conversion 
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Figure D-2.- Concluded 
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Figure D-3*- Planform plot from figure D-2 dataset 
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Figure D-U.- Isometric plot from figure D-2 dataset 
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